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Abstract 


A  higher  order  asymptotic  analysis  of  the  transient  deformation  field  surround¬ 
ing  the  tip  of  a  crack  running  dynamically  along  a  bimaterial  interface  is  presented. 
An  asymptotic  methodology  is  used  to  reduce  the  problem  to  one  of  the  Riemann- 
Hilbert  type.  Its  solution  furnishes  displacement  potentials  which  are  used  to  ex¬ 
plicitly  evaluate  the  near-tip  transient  stress  field.  Crack-tip  fields  corresponding 
to  crack  speeds  up  to  the  loweT  of  the  two  shear  wave  speeds  are  investigated. 
An  experimental  study  of  dynamic  crack  growth  in  PMMA/steel  interfaces  using 
the  optical  method  of  CGS  and  high  speed  photography,  is  also  described.  Tran¬ 
sonic  terminal  speeds  (up  to  1.4c£MMA)  and  initial  accelerations  (~  108m/sec2)  are 
reported  and  discussed.  Transient  effects  are  found  to  be  severe  and  more  im¬ 
portant  than  in  homogeneous  dynamic  fracture.  For  subsonic  crack  growth,  these 
experiments  are  used  to  demonstrate  the  necessity  of  employing  a  fully  transient 
expression  in  the  analysis  of  optical  data  to  accurately  predict  the  complex  dynamic 
stress  intensity  factor  history. 
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1  Introduction 


Advanced  multiphase  materials  such  as  fiber  or  whisker  reinforced  composites  have  seen 
widespread  applications  in  recent  years.  It  has  been  recognized  that  interfacial  fracture 
may  play  an  important  role  in  determining  the  overall  mechanical  response  of  such  multi¬ 
phase  systems.  It  is  the  low  fracture  toughness  of  these  materials,  which  may  result  from 
debonding  between  different  phases,  that  limits  their  use  in  engineering  applications. 
Therefore,  the  scientific  understanding  of  the  mechanics  of  crack  formation,  initiation, 
and  crack  growth  in  bimaterial  interfaces  is  essential  for  the  effective  study  of  failure 
processes  of  these  advanced  composite  materials. 

The  earliest  study  of  interfacial  fracture  appears  to  be  by  WILLIAMS  (1959),  who 
examined  the  local  fields  near  the  tip  of  a  traction  free  semi-infinite  interfacial  crack,  ly¬ 
ing  between  two  otherwise  perfectly  bonded  elastic  halfspaces.  He  observed  that,  unlike 
in  homogeneous  materials,  the  interfacial  crack  exhibits  an  oscillatory  stress  singularity. 
Since  then,  SlH  and  Rice  (1964),  and  RICE  and  SlH  (1965)  have  provided  explicit  expres¬ 
sions  for  the  near-tip  stresses  and  related  them  to  remote  elastic  stress  fields.  The  works 
of  Erdogan  (1965),  England  (1965),  and  Malyshev  and  Salganik  (1965)  have  also 
further  examined  two-dimensional  singular  models  for  single  or  multiple  crack  configura¬ 
tions  in  bimaterial  systems.  Recent  progress  in  static  interfacial  fracture  includes  work 
by  Rice  (1988),  Hutchinson  and  Suo  (1991),  and  Shih  (1991). 

Depending  on  the  nature  of  loads  that  the  composite  structure  is  subjected  to,  the 
debonding  process  may  take  place  dynamically.  If  the  interface  is  already  weakened  by 
pre-existing  flaws,  these  flaws  may  serve  as  sites  of  initiation  of  cracks  which  propagate 
unstably  along  the  interface  under  the  right  circumstances.  Such  situations  have  moti¬ 
vated  attempts  to  analyze  dynamic  crack  propagation  in  interfaces.  However,  due  to  the 
complexity  of  the  problem,  thus  far,  only  a  few  theoretical  results  of  dynamic  bimaterial 
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crack  growth  have  been  obtained.  Among  others,  GoL’DSHTElN  (1967),  Brock  and 
ACHENBACH  (1973),  WlLLIS  (1971,  1973).  and  Atkinson  (1977)  have  provided  crack 
line  solutions  of  particular  fracture  problems.  Although  these  analytical  results  have  pro¬ 
vided  some  insights  of  the  near-tip  dynamic  behavior,  in  order  to  effectively  formulate 
and  apply  crack  initiation  and  growth  criteria  in  bimaterial  systems,  we  need  knowl¬ 
edge  about  the  complete  spatial  structure  of  the  field  surrounding  the  moving  interfacial 
crack-tip. 

More  recently,  experimental  investigations  of  interfacial  crack-tip  deformation  fields 
have  been  carried  out  by  TlPPUR  and  ROSAKIS  (1991)  and  ROSAKIS  et  al.  (1991)  using 
the  optical  method  of  Coherent  Gradient  Sensor  (CGS)  (ROSAKIS,  1993)  and  high  speed 
photography.  The  bimaterial  system  they  used  was  a  PMMA/aluminum  combination. 
They  observed  substantial  crack-tip  speeds  (up  to  90%c^MMA)  associated  with  crack  ini¬ 
tiation  and  growth.  Motivated  by  these  observations,  Yang  et  al.  (1991)  provided  the 
asymptotic  structure  of  the  most  singular  term  of  the  steady  state  elastodynamic  bima¬ 
terial  crack-tip  fields.  In  the  work  of  Wu  (1991),  similar  conclusions  were  reached.  In 
addition,  DENG  (1992)  obtained  the  asymptotic  series  representation  of  the  stress  field 
near  the  tip  of  a  running  interfacial  crack  in  a  bimaterial  system  under  steady  state  con¬ 
ditions.  Also  motivated  by  the  experiments  of  TlPPUR  and  ROSAKIS  (1991),  Lo  et  al. 
(1992)  have  performed  a  numerical  analysis  of  the  same  bimaterial  system  as  was  used 
in  the  experiments. 

The  question  of  whether  there  exists  a  /^-dominant  region  surrounding  the  crack-tip 
(i.e.  a  region  where  the  stress  field  can  be  well  described  by  the  leading  singular  term 
only),  or  in  fact  whether  steady  state  crack  propagation  constitutes  a  good  assumption 
in  analysis,  are  issues  to  be  verified  by  experimental  observations.  New  experimental 
evidence,  described  in  this  paper,  emphasize  the  existence  of  substantial  crack-tip  accel¬ 
erations  in  addition  to  very  high  crack-tip  speeds.  The  existence  of  high  accelerations 
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violates  the  conditions  under  which  the  steady  state  assumption  may  confidently  be  ap¬ 
plied.  Motivated  by  the  above  experimental  evidence,  in  this  paper,  we  investigate  the 
asymptotic  structure  of  the  near-tip  field  in  a  bimaterial  system,  where  a  highly  transient 
elastodynamic  crack  growth  history  has  occured.  To  do  so.  we  employ  the  asymptotic 
procedure  proposed  by  FREUND  (1990)  and  utilized  by  FREUND  and  ROSAKIS  (1992)  in 
studying  the  transient  growth  of  a  mode-I  crack  in  a  homogeneous  isotropic  material  The 
same  procedure  was  employed  by  LlU  and  ROSAKIS  (1992)  in  studying  the  mixed-mode 
transient  growth  of  a  crack  along  an  arbitrary  curved  path  in  a  homogeneous  isotropic 
solid.  For  anisotropic  solids,  transient  crack  growth  under  mode-I  conditions  was  recently 
explored  by  WlLLIS  (1992). 

In  section  2  of  the  present  study,  the  general  formulation  and  properties  of  the  asymp¬ 
totic  procedure  are  described.  By  using  this  asymptotic  methodology,  the  equation  of 
motion  is  reduced  to  a  series  of  coupled  partial  differential  equations.  In  section  3,  the 
solution  for  the  higher  order  transient  problem  is  obtained.  By  imposing  the  boundary 
conditions  along  the  surface  of  the  interfacial  crack  and  the  bonding  conditions  along  the 
interface  ahead  of  the  crack-tip,  our  problem  can  be  further  recasted  into  a  Riemann- 
Hilbert  problem.  Upon  solving  the  Riemann-Hilbert  equation  and  evaluating  the  Stieltjes 
transforms,  the  higher  order  near-tip  transient  elastodynamic  asymptotic  field  can  be  ob¬ 
tained.  In  section  4,  the  asymptotic  elastodynamic  stress  field  surrounding  the  interfacial 
crack-tip  is  studied.  The  first  stress  invariant  is  provided  explicitly.  The  properties  of 
the  interfacial  mismatch  parameters  are  studied  in  section  5.  These  depend  on  the  prop¬ 
erties  of  the  bimaterial  combination  and  the  crack-tip  speed.  In  some  of  the  available 
experiments  by  ROSAKIS  et  al.  (1991),  and  the  experimental  evidence  described  in  this 
paper,  it  has  been  observed  that  an  interfacial  crack  can  reach  speeds  amounting  to  a 
considerable  fraction,  or  even  exceeding  the  lower  Rayleigh  wave  speed  of  the  two  con¬ 
stituents  of  the  interface.  Recognizing  that  our  analysis  need  not  be  limited  to  a  velocity 
regime  below  the  lower  Rayleigh  wave  speed,  in  section  6,  we  extend  our  solution  to  the 
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case  when  the  crack  is  travelling  at  a  speed  between  the  lower  Rayleigh  and  shear  wave 
speeds.  Finally,  in  section  7,  recent  experimental  evidence  of  a  transient  higher  order 
stress  field  in  bimaterial  fracture  specimens  is  presented.  The  transient  theoretical  fields 
obtained  in  previous  sections  are  used  to  quantitatively  analyze  optical  interferograms 
obtained  in  real  time  high  speed  photography  of  dynamic  birnateria!  experiments  in  a 
PMMA/steel  system.  In  addition,  we  present  experimental  evidence  of  transonic  crack 
growth  histories  involving  maximum  speeds  between  60%  and  80%  of  the  dilatational 
wave  speed  of  PMMA.  For  comparison  purposes,  one  should  note  that  typical  terminal 
crack-tip  speeds  in  homogeneous  PMMA  are  of  the  order  of  only  20%  of  the  dilatational 
wave  speed. 

2  General  formulation 

Consider  a  planar  body  composed  of  two  homogeneous,  isotropic,  and  linearly  elastic 
materials  which  are  bonded  along  a  straight  line  interface.  A  crack  propagates  non- 
uniformly  along  the  interface,  see  Figure  1.  Introduce  a  fixed  orthonormal  Cartesian 
coordinate  system  (xi,x2)  so  that  the  xj-axis  lies  on  the  interface  and  coincides  with 
the  direction  of  the  propagating  crack.  Suppose  that  the  crack  propagates  with  a  non- 
uniform  speed,  u(<),  and  the  crack  faces  satisfy  traction  free  boundary  conditions.  At  a 
time  t  =  0,  the  crack-tip  happens  to  be  at  the  origin  of  the  system,  so  the  growth  of  the 
interfacial  crack  at  any  t  >  0  is  characterized  by  the  length  l(t)  (v(t)  —  which  is 

the  distance  from  the  origin  to  the  moving  crack-tip.  If  the  deformation  is  assumed  to 
be  plane  strain,  for  each  of  the  two  materials  comprising  the  interface,  the  displacement 
field  may  be  generated  from  two  displacement  potentials,  <f>k(xux2,t)  and  0fc(xi,x2,t), 
where  k  €  {1,2}.  Here,  the  integer  k  is  assigned  to  distinguish  between  the  two  different 
materials.  In  Figure  1,  material-1  is  the  one  shown  above  while  material-2  is  the  one 
shown  below  the  interface.  Then,  in  either  one  of  the  two  materials,  the  two  non-zero 
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displacement  components  can  be  expressed  as 


,  -^2,  0  —  0,a  (^1 »  % 2,  t )  +  Cai 30, 0  (-^l ,  %2i  ,  (2. 1  ) 

where  a,  fi  €  {1,2}  and  the  summation  convention  has  been  used.  ea$  is  the  two- 
dimensional  alternator  defined  by 


e12  —  ~  ^21  —  1  ,  en  —  ^22  —  0  • 


The  components  of  stress  for  each  material  can  be  expressed  in  terms  of  the  displacement 
potentials  by 


On  =  n 


c2 

“70, ao,  ~ 20,22  +20, 12 


(P" 

L  S 


^22  —  A4  “70, era  ~20,1i  —20,12 


(2-2) 


°12  —  A4  |^20,12  +0,22  ~ 0,11 1  J 

where  a4  is  the  shear  modulus,  and  c/,  c,  are  the  longitudinal  and  shear  wave  speeds  of 
each  elastic  material  above  or  below  the  interface,  respectively.  In  terms  of  the  shear 
modulus  a4,  nass  density  p,  and  Poisson’s  ratio  i/,  for  each  of  the  two  materials,  c/  and 
c,  are  given  by 

where 

3  —  4i/,  Plane  strain 
*  =  ^  3-«/ 


1  +i/’ 


Plane  stress 


The  corresponding  plane  stress  solution  can  be  obtained  by  changing  the  definition  for 
the  longitudinal  wave  speed  in  equation  (2.3).  Meanwhile,  q  and  c,  in  both  plane  strain 
and  plane  stress,  are  related  by 


Cm  _  (K_ 

q  Ik 


c,  r  *  - 1 1 1/2 
+  i  / 


(2.4) 
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The  equation  of  motion  in  the  absence  of  body  forces  in  the  fixed  coordinate  system, 
in  terms  of  <^(xj,x2,f)  and  ^>(xi,x2,<),  reduces  to 


<t>,aa{xux2,t) - r^(xi,x2,t)  =  0 

ci 

Vw(xi,x3,<)  -  -^0(xi,x2,<)  =  0 


> 


(2.5) 


Equations  (2.5)  hold  for  each  material  above  or  below  the  interface. 


We  further  introduce  a  new  moving  coordinate  system,  (£ j  ,  £2),  by 


6  =  ~  J(<)  , 


6  =  *2  • 


(2.6) 


This  system  is  such  that  its  origin  is  translating  with  the  crack-tip.  In  this  new  system, 
the  equations  (2.5)  for  0(6,6,  t)  and  0(6,6,*)  become  (Freund,  1990) 


(-?) 

(-?)*« 


0,11  +0,22  +  ““5^0,1  H - 4"^0,1<  -~<f>,tt  =  0 


C? 


-2 

c/ 


L-; 


,  t>(<)  ,  2u(<)  .  1  .  . 

+0,22  H - “0,1  H - - 0,11 - j0,tt  =  ® 

c,  c.  c. 


(2-V) 


Notice  that  in  equations  (2.7),  the  differentiation  with  respect  to  time,  t,  is  distinct  to  that 
in  equations  (2.5).  Here,  (6,6)  are  held  fixed,  whereas  in  (2.5),  (ii,x2)  are  held  fixed. 
Throughout  this  study,  we  will  use  dfdt ,  or  {  },<  to  denote  differentiation  with  respect  to 
time,  t,  when  the  moving  spatial  coordinates  (6,6)  are  held  fixed.  The  notation  {  } 

denotes  the  same  operation  when  the  spatial  coordinates  (xj,x2)  are  held  fixed. 

At  this  point,  we  employ  the  standard  asymptotic  device  used  by  Freund  and 
ROSAKIS  (1992)  for  the  analysis  of  transient  crack  growth  in  homogeneous  materials. 
We  assume  that  0(6,6, 0  and  0(6,6,*),  for  each  material,  can  be  asymptotically  ex¬ 
panded  as 

OO 

0(6,6. 0  =  £ePm0m{iJ  1,72,0 


m=0 

OO 


0(6,6, 0  =  HePm0m('?i 


m=0 


,9a,*)  J 


(2.8) 
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as  r  =  (£1  4-  £f)^2  -+  0,  where  t)q  =  £a/e  ,q  €  {1,2),  and  £  is  a  small  arbitrary  positive 
number.  The  parameter  £  is  used  l.ere  so  that  the  region  around  the  crack-tip  is  expanded 
to  fill  the  entire  field  of  observation.  As  &  is  chosen  to  be  infinitely  small,  all  points  in 
the  (£1,^2)  plane  except  those  very  close  to  the  crack-tip,  are  pushed  out  of  the  field 
of  observation  in  the  (771,772)  plane,  and  the  crack  line  occupies  the  whole  negative  t?i - 
axis.  By  taking  e  =  1,  the  above  equations  will  provide  the  asymptotic  representation 
of  the  displacement  potentials  in  the  unsealed  physical  plane  for  each  of  the  materials, 
respectively. 


In  the  asymptotic  representation  (2.8),  the  powers  of  s  are  such  that 

Pm+l=Pm  +  ^,  771=0,1,2,  •••.  (2.9) 

Since  the  displacement  should  be  bounded  throughout  the  region,  but  the  stress  may  be 
singular  at  the  crack-tip,  p0  is  expected  to  be  in  the  range  1  <  po  <  2  (FREUND,  1990). 
We  also  should  have 


£pm<t>m(v  1,W) 


(2.10) 


for  any  positive  integer  n.  Returning  to  the  unsealed  physical  plane,  we  will  have 


^m4n(£l,  £2,0 

^m(£l,£2,  0 


as  r  -  \f[\  +7f  -»  0  , 


(2.11) 


for  any  positive  integer  n,  so  that  in  the  physical  plane,  (£i,£2),  <f>m(£i,£2,0  are  ordered 
according  to  their  contributions  to  the  near-tip  deformation  field.  The  above  properties 
for  4>m  hold  for  0m  as  well. 


Substituting  the  asymptotic  representations  for  <tK£i>£2>*)  an<^  0(£ i>£2,0>  (2-8),  into 
the  equations  of  motion,  (2.7),  we  obtain  two  equations  whose  left  hand  side  is  an  infinite 
power  series  in  £  and  whose  right  hand  side  vanishes.  Since  e  is  an  arbitrary  number,  the 
coefficient  of  each  power  of  e  should  be  zero.  Therefore,  the  equations  cf  motion  reduce 
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to  a  series  of  coupled  differential  equations  for  <pm{r)  1,172,  t)  and  ,  J?2-  0  as  follows 

(Rosakis  et  a/.,  1991,  FREUND  and  ROSAKIS,  1992): 


1  1  1 
^m,n  +^f(7)^Tn'22  =  ”  af(0c?  K  +  I 

0m’n+^7)^’22  =  “Snw  l^^-2’1  K + 


(2.12) 


»2(0  “  a2(<)C2  1"  v','rr"'‘,‘  i'i  '  a2(t)( 

for  m  =  0, 1,2,  •••,  and  the  quantities  q;  and  a,  depend  on  the  crack-tip  speed,  and 
therefore  on  time  t  through 


°l(i)  =  1  - 


v2(t) 


-1,, 


Also 


4>n 


■{ 


4>m  for  m  >  0 
0  for  m  <  0 


— 


I 


i/>m  for  m  >  0 
0  for  m  <  0 


(2.13) 


(2.14) 


In  what  will  follow,  for  our  convenience,  we  drop  the  subscript  which  is  used  to 
distinguish  between  the  two  materials.  However,  we  should  keep  in  mind  that  the  above 
asymptotic  form  of  the  equations  of  motion  (2.12)  hold  for  each  of  the  materials  with 
the  appropriate  elastic  constants.  The  term  “coupled”  is  used  above  in  the  sense  that 
the  higher  order  solutions  for  <j>m  and  ipm  will  depend  on  the  lower  order  solutions  for 
the  same  quantities.  It  is  noted  that,  for  the  special  case  of  steady  state  crack  growth, 
the  crack-tip  speed,  t>,  will  be  a  constant,  and  at  the  same  time,  Vw  =  0,  for 

m  —  0, 1,2,  •  •  •.  This  means  that  4>m  and  ipm  depend  on  t  only  through  the  spatial  scaled 
coordinate  ifo.  In  such  a  case,  the  equations  in  (2.12)  are  not  coupled  anymore,  and  each 
one  reduces  to  Laplace’s  equation  in  the  coordinates  (71,07772)  for  <j>m  and  (7 1,0,72)  for 
tpm.  For  steady  state  conditions,  the  functions  <f>m  and  are  independent  of  time  in 
the  moving  coordinate  system.  For  the  transient  case,  however,  the  crack-tip  speed,  u(f), 
may  be  an  arbitrary  smooth  function  of  time,  and  also  4>m  and  ipm  may  depend  or  time 
explicitly  in  the  moving  coordinate  system.  The  only  uncoupled  equations  are  those  for 
m  =  0  and  m  =  1.  For  m  >  1,  we  can  see  from  (2.12)  that  the  solutions  for  <£m,  and  t/>m 
are  composed  of  two  parts.  One  is  the  particular  solution  which  is  whoiely  determined 
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by  lover  order  solutions  for  <f>,n,  and  The  other  part  is  the  homogeneous  solution 
which  satisfies  Laplace’s  equation  in  the  corresponding  scaled  coordinate  plane.  The 
combination  of  the  particular  and  homogeneous  solutions  should  satisfy  the  traction  free 
conditions  on  the  crack  faces  as  well  as  the  bonding  conditions  along  the  interface,  in  the 
following  sections,  we  will  solve  for  <pm  and  4>m  for  the  most  general  transient  situation, 
and  for  both  materials. 

It  should  be  noted  that  the  steady  state  problem  could  be  solved  using  the  efficient 
Stroh  formulation.  This  formulation  reduces  the  two  spatial  and  one  temporal  variables 
to  only  two  spatial  variables  and  takes  advantage  of  a  well  known  formalism  to  solve  the 
steady  state  crack  problem  (Yang  et  al.,  1991).  However,  this  approach,  although  it  can 
be  easily  be  extended  to  anisotropic  solids,  is  strictly  restricted  to  steady  state  conditions 
and  cannot  be  used  for  our  present  purposes. 
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3  Solution  for  the  higher  order  transient  problem 


As  we  have  discussed  in  the  previous  section,  the  only  uncoupled  equations  in  (2.12)  are 
those  for  m  =  0  and  m  =  1.  For  m  >  1,  the  solutions  for  <t>m  and  will  be  affected 
by  the  solutions  with  smaller  m.  In  this  section,  we  consider  the  situation  of  m  =  0  and 
m  —  1  first.  After  we  get  solutions  for  m  =  0,  and  1,  we  will  subsequently  solve  for  higher 
order  4>m  and  t/>m. 


3.1  Solutions  for  *?2,  t)  and  372, t)  for  m  =  0  and  1 


For  m  =  0,  or  1,  (2.12)  reduce  to 


0m. 11  (»?li»?2i0  d - 277T<^m,22  (*7l  )  *?2.  0  =  0 

Q,(<) 

<*;(*) 


(3.1) 


The  above  equations  are  Lapalace’s  equations  hi  the  corresponding  scaled  planes  (*h,  cqr?2) 
for  <f>m,  and  (»7i, 0,172)  for  ifim.  As  we  have  mentioned  earlier,  the  subscript  k  is  omitted 
here,  but  the  above  equations  hold  for  both  materials  that  constitute  the  bimaterial  body. 


The  most  general  solutions  for  equations  (3.1)  can  be  expressed  as 

0m(*?i,»?2,O  =  Re{Fm(z,;f)}  ,  =  Im{Gm(z,;t)}  ,  (3.2) 

where  the  two  complex  variables  zi  and  z,  are  given  by 

*1  =  m  +  *ari»fe  ,  z$  =  rji+  iatrj2  , 

and  i  =  v—1.  For  the  bimaterial  system,  Fmk(zr,t)  and  Gmk(z,;t)  are  analytic  in  the 
upper  half  complex  z/*- ,  or  z,*-planes  for  k  =  1  (upper  material),  and  analytic  in  the 
lower  half  complex  zj*-,  or  z4*-planes  for  k  —  2  (lower  material).  The  complex  conjugates 
of  these  functions  are  also  analytic  on  the  plane  obtained  by  reflection  along  the  real  axis, 
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e.g.  Fmi(li\t)  is  an  analytic  function  on  the  zt  plane.  Since  as  and  a,  are  different  for 
each  material,  the  scaled  complex  variables  z\  and  z,  will  also  be  different.  For  fully 
transient  problems,  in  the  analytic  functions  Fmk(zt\t)  and  Gmk{za\  t ),  time  t  appears  as 
a  parameter.  This  suggests  that  these  functions  will  depend  on  time  t  not  only  through 
the  complex  variables,  zi  and  zt,  but  also  directly  through  time  t  itself. 


The  displacement  and  stress  components  associated  with  these  <pm  and  rpm,  are  given 


u(r)  =  Re{^(,i;0  +  aJG;(z,;0} 

(3-3) 

up1  =  . 

and 

<*’  =  ,,Re{(l  +  2a?  -  a?)  F^z,-,  t )  +  ()}  ' 

a<r>  =  -A.Re{(l+a;)F"(2,;i)  +  2a.G^(2.;()}  ;  .  (3.4) 

=  -pirn  {2a, +  (l  +  a?)  C"  (*,;!)} 

where  primes  denote  the  derivative  with  respect  to  the  corresponding  complex  arguments. 


For  any  analytic  function  fi(z),  define  the  following, 

lim  fl(z)  =  fl+(»7i) 
o+ 

’  ,  ^  m  • 

lim  fl(z)  =  n~(T?i) 
m— 0“ 

For  t)i  <  0  and  t)2  — *  0+,  the  traction  free  condition  on  the  upper  crack  face  gives 

{<£’},  =  {dr1},  =  0  , 

or,  in  terms  of  the  complex  displacement  potentials  Fm(zi;t)  and  G„,(zf;<), 

{/>(1  +  «J)  [C+(„;<)  +K („;<)]  +  2?a.  [g«(,i;()  +  G,:-(„;()]}1  =  0 

I  ’  (3’5) 

{2 nm  -  C On; <)]  +  /»(i  +  «.!)  [<C0n; «)  - BT(i.; <)] },  =  o 

For  r?i  <  0  and  the  traction  free  condition  on  the  lower  crack  face  gives 

{<£’},  =  {»!?’}, = o , 


II 


or 


{Mi  +  oj)  [C (■?.;<)  +  Oli; <)]  +  2/iq,  0  +  On; <)] },  =  o  | 

}  (3.6) 

{2/10;  -  On;')]  +M1  +  °1)  |<C (n:»)  -On; <)]},  =  0  j 

The  above  equations,  (3.5)  and  (3.6),  hold  for  rji  <  0. 


Along  the  interface,  or  ifo  >  0  and  =  0,  the  bonding  conditions  should  be  satisfied, 
which  implies  that 


Mr*}, -{<£'}, -o,  {<#>},  -  {<#’},  =  o 

RVKV*  H-’h-m.-o, 


V  rji  >0.  t)2  —  0  , 


or,  in  terms  of  Fm{zi\t)  and  Gm(z,]t),  from  traction  continuity, 

{Ml  +  «;)[On;0  +  On;<)]  +2/<a,  [G7(n;0  +  On;<)]}, 

-{Mi  +<*:)  [C  (n;<) +On;0]  +2/-0.  [On/O+On;')]},  =  o 

, 

{*m  \K+(m-,t)  -TT(«h;0]  -Md  +  «!)  [^(wO-iCtoi;*)]}, 

-{2jra,  [F"'(j7!;t)  -  F^+(j?,;0]  +  M1  +<*2)  “^+(»?i;0]}a  =  0  , 

(3.7) 


and  from  displacement  continuity, 

{ [On;  0  +  K  (n ;  <)]  +  o,  [On  it) + On ;  o] }, 

-  {[C(>h;<)  +  F^(7.;<)l+o.[C(’Jii<)  +  G'n;('ii;<)]}2  =  o 

{o,  [On;<)-On;<)]  +  [On;0-On  ;<)]}, 

-  {a,  [On;0-On;<)]  +  [On;') -On;')]},  =  o  , 

The  above  equations,  (3.7)  and  (3.8),  hold  for  i?i  >  0. 


(3.8) 


For  simplicity,  define  the  following  matrices  for  each  material,  k  £  {1,2}, 


Pk 


n{\  -I-  a])  2 fiat 

2fiat  n(l  +  a*)  J 


Qk  - 


/i(i  +qJ) 

-  2pOq 


2  fiat, 

-Ml  +  <*?)  J 
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uk  = 

1 

Q.  " 

II 

■46 

1  Q» 

.  Ql 

1 

k  ! 

-  a,  -l 

Also  define  the  following  complex  vector  for  each  material, 

Gmk{z;t)  )T  , 


where  z  =  rj i  +  ir/2-  From  above  definitions,  the  boundary  and  bonding  conditions, 
equations  (3.5),  (3.6),  (3.7),  and  (3.8),  can  be  rewritten  as 


Pifmiin  i;  0  +  Qj'miiv  1;  0  =  ° 


v  »?i  <  o  , 


(3-9) 


and 


Pif^vut)  +  QJT, (*;«) - =  o  ] 

1 ,  V  i),  >  o.  (3.10) 

t'.Oi.;')  +  v'.T'J.to.K)  -  u.f-^t)  -  =  o  ) 


Further,  the  bonding  conditions  in  (3.10)  can  be  rearranged  as 

-  Qjfcfat)  =  *»/£(*; 0  -  QiTHTiOn;*)  ' 

*  > 

Uif'+Mut) -  v27+,(n;t)  =  -  v,7^(m-,t)  . 


Viji  >  0.  (3.11) 


In  the  above  equations  (3.11),  the  left  hand  sides  are  the  limiting  values  of  functions 
which  are  analytic  in  the  upper  half  plane.  The  right  hand  sides  are  the  limiting  values 
of  functions  which  are  analytic  in  the  lower  half  plane.  Since  the  limiting  values  are  the 
same  along  the  positive  real  axis,  the  function  Pifmii2^)  ~  Qifmuzi*)  defined  in  the 
lower  half  plane,  is  the  analytic  continuation  of  the  function  Pif^z^t)  — 
which  is  defined  in  the  upper  half  plane,  and  vice  versa.  This  results  from  the  continuation 
propertities  of  analytic  functions.  As  a  result,  we  cam  write 


^ 0  —  Ql/ml(z>  0  =  Km(zi  0  7  2  6  S 


(3.12) 
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and 


2  €  S+ 


(3.13) 


Uif^frt) -  vj'ml(z;t)  =  em(z;t) ,  zes- 


where 


S* 


|  {  (Vu  m)  1  -  oo  <  Tjy  <  OO,  T)2  >  0  }  -  C 


{  {  (vu  ”72)  |  -  00  <  J?1  <  OO,  T)2  <  0  }  -  C 

c  =  {  (»?1,  irj2)  j  -  00  <  J?!  <  0,  1)2  =  0  }  . 

Km(z\ f)  and  8m(z',t)  are  analytic  functions  throughout  the  2-plane  except  along  the  cut 

C  which  is  the  entire  non-positive  real  axis.  From  the  above  equations,  it  can  be  seen 

immediately  that  equations  (3.10)  are  satisfied  identically.  So,  the  question  now  is  to  find 

the  analytic  functions  Km(z;t)  and  Om(z;t)  in  the  cut-plane  S+  U  S~ . 


Solving  for  f'^k(z‘,t)  and  f'^k(z;t)  from  equations  (3.12)  and  (3.13),  we  get 

=  er'ff-1  {•'.(*;<)-  U  *-(*;<>} 


2€  5+  , 


(3.14) 


and 


/:,(*;<)  =  -PV  H  i,  «»(*i0} 

71i(*;  t)  =  -Or1  h'  t )  -  ()} 


2  6  5’”. 


(3.15) 


The  definitions  of  matrices  £*,  Lk ,  H,  and  H,  as  well  as  the  properties  of  these  matrices 
are  given  in  Appendix  I.  Matrices  Pk  and  Qk  have  been  defined  above.  In  obtaining 
(3.14)  and  (3.15),  we  have  assumed  that  the  inverse  matrices  Pkl  and  Qk 1  exist.  Notice 
that  the  determinants  of  Pk  and  Qk  are  both  equal  to  Dk(v),  where 

£»i(v)={4a/Q,-(l  +  a,J)2}t  . 

Therefore,  in  this  analysis,  we  exclude  the  situation  where  the  interfacial  crack  propagates 
with  either  of  the  two  Rayleigh  wave  speeds  of  the  bimaterial  system  which  are  the  real 
roots  of  Dk(v)  =0.  This  ensures  the  existence  of  Pk  1  and  Qkl. 


14 


Substituting  equations  (3.14)  and  (3.15)  into  the  traction  free  conditions  on  the  crack 
faces,  (3.9),  we  get 

H  ^2  -  H  {OwO  -  LlK-miVut)}  =  O 

<0.  (3.16) 

H  t)-  Lx  K-  (*;<)}  -  H  {•£(*»;  t)  -  £,<(*,;  *)}  =  o  j 

m 

Adding  the  two  equations  in  (3.16),  and  using  the  fact  that  H  H  /  o  for  a  crack 
propagating  with  a  nonzero  speed,  we  obtain 

«£0h;0-*m0?i;0  =  ° »  Vr?i<o.  (3.17) 

This  implies  that  Km(z\t)  is  continuous  across  the  cut  except  at  the  crack-tip  and  there¬ 
fore  Km(z;t )  is  analytic  in  the  entire  complex  plane  except  at  z  =  0.  However,  the  con¬ 
dition  of  bounded  displacement  requires  that  |  k m(z;  f)  |  =  O  (  \z\a  )  for  some  a  >  —  1, 
as  \z\  — *■  0,  so  that  any  singularity  of  Km(z;t)  at  the  crack-tip  is  removable.  Therefore, 
Km(z;t)  is  an  entire  function.  Now,  both  equations  in  (3.16)  become 

H  ,  V,1<0.  (3.18) 

where 

*m(*n;0  =  <(ii  -,t)  =  *m(v  i;0 

R  =HL2  -HL2  =H  Li  -H  i, 

Equation  (3.18)  constitutes  a  Riemann-Hilbert  problem.  Its  solution  B'm{z\t)  is  ana¬ 
lytic  in  the  cut-plane  S+  U  S~.  Along  the  cut,  0'm(z-,  t)  satisfies  equation  (3.18)  for  some 
arbitrary  entire  function  Km(z;t).  Also,  from  the  requirement  of  bounded  displacements 
at  the  crack- tip,  as  |z|  — ►  0, 

I4.(*;0l-O(l*n,  (3.19) 


for  some  a  >  — 1. 


Id  equation  (3.18),  the  solution  6'm(z\t)  is  composed  of  two  parts,  the  homogeneous 

-  t 

solution  6m  ( z;t ),  and  the  particular  solution  0m(z]t).  We  will  consider  these  two  solu¬ 
tions  separately. 

Homogeneous  Solution: 

The  homogeneous  solution  is  obtained  by  solving 

(t?i;  t)  -  H  0m  (rjut)  =  o  ,  V  ^  <  0  .  (3.20) 

By  using  the  solution  given  in  Appendix  II  and  by  imposing  restriction  (3.19),  we  can 
write  the  solution  to  the  above  equation  as  follows: 

K  (*;<)  =  2-U"  °Am  (*;  !)C  +  Bm  (z;  t)  C  ,  (3.21 ) 

where 

2*1  +  /?’  P~  y/knh^' 

<  =  (i,  »?)T>  C=(i, r)  =  )[^2 

o  ° 

and  Am  (z;<)>  J&n  (2;  t)  are  arbitrary  entire  functions.  The  parameters  e  and  tj  defined 
here  are  known  functions  of  crack-tip  speed,  v,  and  material  properties.  Their  functional 
dependence  on  these  variables  is  discussed  in  Appendix  I  and  section  5.  For  v  =  0,  e(u) 
reduces  to  Cq  which  is  the  oscillatory  index  that  appears  in  the  quasi-static  interfacial 
crack  problems  (WILLIAMS,  1959;  Rice,  1988). 

By  substituting  equation  (3.21)  into  equations  (3.14)  and  (3.15),  we  get 

>1,  (*; 0  =  j Vtf-1  {z“i+“  °Am  (z:() f  +  z-i-”  iu  <z; f)  <}  | 

>  ,  2  €  5+ ,  (3.22) 

/„,(*;<)  =  Q^H-’  {z*U“  Am  (z;()C  +  Z-U”  e„  (z;t)  c}  j 
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and 


)m2  (z;t)  =  -P,1  H  {z-$+“Am(z;tK  +  z-t-'(  Bm(z;t)c}  j 

/mt(*;  t )  =  -Qr1  H1  {*-*+•■'  Am  (z;  t)c  +  Bm  (2;  0  c}  J 


},  z€5“.  (3.23) 


Notice  that  the  following  identities  hold, 

\Ai2Ml 


Ml  —  ^12^21 


(i-/J)f,  h-1  C=  Tr^T^f-(i  +  fi)  < 


a"  c  =  ir^rVt1  +  /»)C,  «"c= 

«ll  ~  «12«21 


Ml  -  M2M1 

12^21 

Ml  —  M2M1 


(1-flC 


d- 


and 


l+/3  = 


1-/? 


cosh  «r  ’  cosh  tv 

Without  losing  generality,  we  may  absorb  the  factor  \Ai2Mi/(Mi  -M2M1)  into  the  entire 
00  T" 

functions,  Am  (2;*)  and  J3m  (z;t).  By  taking  the  conjugate  of  the  function  /ml(z;t)  in 

o" 

equation  (3.23)  and  comparing  it  with  the  function  /ml  ( z ;  <)  in  equation  (3.22),  and  also 
by  using  the  properties  of  matrices  Pk  and  Qk,  we  can  obtain  a  relationship  between  the 
entire  functions  Am  (2;  t)  and  Bm  (2;  t)  as  follows, 


Bm  (2;  #)  =  -Am(z;t)  ■ 


Meanwhile,  by  using  the  fact  that 


p;'  c  =  q;'c  ,  <?;’<=  p;'c  , 


we  can  get  the  solutions, 

o»  f  e**P~xC  1  o 

(*i  *)  =  — £-M+"  A.  (*;<)  + 


cosh  eir 


Am  (2 


€  S+ 


°fm2  cosh£2rCz'i+,t^m(z;<)+  co!h2e7rC2'i'U^m(z;f)}’  Z€5~  , 


(3.24) 
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or,  in  terms  of  Fm(zr,t)  and  Gm(z,\t),  for  ihe  material  above  the  interface, 


Gm(z,-,t)  = 


[(l  +  aJ2)-2r?aa]et'  _i+,{  •  , . 

- ^0(i>)cosheir .  *' 

1(1  +al)  +  27}aa}e-‘* 

- „!>(„)  cosT”'2'  Am(2ht) 

[2Q/~t?(  1  +  a,)]  e<ir  ~^+,t  ° 

nD(v)  cosher  *  m  *’ 

[2a/  +  r?(l  +  a*)]e 


—tir  |  .  rr 

-4-u  0 


+  —  nr;  —r1 — 1 

fiD(v)  cosh  eff 


(3.25) 


For  the  material  below  the  interface,  the  solution  is  also  given  by  equation  (3.25)  with 
the  parameter  e?r  changed  to  —en. 


Particular  Solution: 

Since  *cm(z;t)  is  an  entire  function,  the  particular  solution  0m(z;t)  can  be  easily 

A  /  ^  t 

constructed.  Suppose  Om(z;t)  is  also  an  entire  function,  which  implies  that 

then  from  equation  (3.18),  we  get 

«!(>?.;<)=  {A-h}’1  RK„(Vl;t),  V„<0.  (3.26) 

By  using  the  identity  theorem  for  analytical  functions,  it  can  be  shown  that  for  any  z, 

Oliz-J)  =  {k -hY' RKm(z;t)  .  (3.27) 


By  substituting  this  particular  solution  into  equations  (3.14)  and  (3.15),  we  have 


/!,(*; <)  =  p;'  {h  -h}“  {i2  -£,}  t) 

7U’-,t)  =  -Qi'  {h  (t,  -L,} 


Z  €  S+ , 


(3.28) 
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and 


?M(z-,t)  =  p;' [h -h}'  (i, 

%, (*;«)  =  -QT'  |/r  {£,  -£,}  «_(*;«) 


2  €  S' 


Notice  that 


(^21  )t 


{h-h}"  {£*-£*}=  hn 


(/«)*  ’ 


fc  €  {1,2} 


(3.29) 


If  the  entire  function  /cm(2;  t)  is  expressed  as 


*-.(*;*)  =  (  *£}(*;0  )T  , 

then,  it  can  be  shown  that  by  comparing  the  conjugate  of  fmX{z’,t)  in  equation  (3.29) 

*  ft 

with  in  equation  (3.28),  we  have 

0  +  ^m(zi  0  =  0,  0  -  0  =  0. 


Define  a  new  entire  function  Am{z\t)  by 

Am(*;0  =  t{[«»  (*;0  -  «m(2i0]  +  [*m)(**»0  +  *52)(*i0]} 


Also  let 


0ai)fc  ihi)k 

h2i  h\2 


By  relating  Km(z;t)  to  v4m(z;t),  and  by  using  the  above  definition,  equations  (3.28)  and 


(3.29)  give 


f"mX{z-J)  =  P;'w2Am(z-,t)  -  Q~lw2Am(z-,t)  ,  z  €  5+ 
f"m 2(2; 0  =  P?tvxAm{z\t)  -  QilwxAm(z;t )  ,  2  €  S' 


(3.30) 
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In  order  to  express  the  particular  solution  in  terms  of  Fm(zi\  t)  and  Gm(z,;  t),  we  need 
to  define  two  parameters,  to,  and  to,  that  only  depend  on  the  crack-tip  speed, 


to, 


fa,(l  -  of)' 

fa,(l  -  of)' 

l  vD{v)  j 

i  1  t*D(v)  . 

L 

[a,(l  of) ' 

fo,(l-of)l 

1  »D(v)  . 

(;  1  vD(v)  J 

t 

Then,  for  the  material  above  the  interface,  the  particular  solution  can  be  expressed  as 


<&(*.;<)  = 


1 


\  (  2a,  1  +  of  , 

|(j77r  -TT“ 


fiD(v)  ^  1^1  -f  to,  1  +io, 


(3.31) 


For  the  material  below  the  interface,  the  particular  solution  is  also  given  by  equations 
(3.31)  with  to,  and  to,  changed  to  tof1  and  to”1,  respectively. 

By  adding  the  expressions  in  equations  (3.25)  and  (3.31),  and  by  integrating  with 
respect  to  the  corresponding  arguments,  the  final  solutions  of  Fm(zt;  t )  and  Gm(z,;  t )  for 
the  material  above  the  interface,  for  m  =  0,  1,  can  be  obtained  as 


Fm(z,;t) 


1(1 -fa,2)- 2*70,] 


p/>(u)coshcjr 


zf  Am{zr,t) 


[(1  +  of)  +  2*70,]  e~ 


a_ 


zf  Am{zr,t ) 


PD{ 


pD(v)  cosh  c*r 


(3.32) 
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and 


n  /_  .  _  [2q;  *?(1  +  Qa)]g  f+lf4  /  . 

Gm(  *’  )  ”  MD(v)  cashew  *  AmM) 

■  [2Qj  +  7y(l+Q;)]e-^ 

+  fiD(v)  cosh  C7T  2*  m(  ,,  ) 


(3.33) 


where  the  entire  functions,  Am(z\  t)  and  Bm(z;t )  are  defined  by 


j^{z^Am(z;t)}  =  z~±+i‘  Am  (z;t)  ,  ~{z2£Uz;f)}  =  , 


and  they  can  only  be  determined  by  the  far  field  conditions.  The  solutions  for  the  two 
displacement  potentials,  <t>m(7h,V2,t)  and  tj2,  <),  will  be  given  by  equations  (3.2). 


Since  Am{z\t)  and  Bm(z;t)  are  entire  functions,  they  can  be  expanded  into  Taylor 
series, 

m*  t)  =  £;  4nV)*n .  b0[z; t) = f; 

n=0  n=0 

M*;  <)  =  L  =  f;  B'-’Mz" 

n=0  n=0 

As  we  have  mentioned  in  the  previous  section,  in  the  unsealed  physical  plane,  (£1,(2), 
and  ipm(( j,£2»0  should  be  ordered  according  to  their  contributions  to  the 
near- tip  deformation  field.  By  imposing  this  property,  i.e.  equation  (2.11),  to  the  repre¬ 
sentations  of  ^m(*h and  for  m  =  0  and  1,  we  can  obtain  restrictions 

on  the  entire  functions  Am(z;  t)  and  Bm(z-,t).  In  the  Taylor  expansion  (3.34),  Ao0)(t)  ^  0 
and  i?{0^(f)  7^  0,  but  A^(t)  =  0.  In  other  words,  the  leading  terms  of  fa  and  tl>0  are 
of  order  z3/2,  whereas  the  leading  terms  of  <f>  1  and  $>\  are  of  order  z2.  Meanwhile,  it 
can  be  shown  that  the  coefficient  of  the  leading  term,  Ab°\t),  in  (3.32)  and  (3.33),  is 
directly  related  to  the  complex  dynamic  stress  intensity  factor  Kd(t)  defined  by  Yang  et 
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al.  (1991)  through  the  relation 


4o)(0  = 


_1 _ K\t) 

^8 -HI*-)' 


(3.35) 


As  a  matter  of  fact,  in  the  unsealed  plane,  ((i,(2),  and  for  m  =  0,  equations  (3.32) 
and  (3.33)  are  identical  in  spatial  structure  to  the  complete  solution  for  the  steady  state 
propagating  interfacial  crack  in  a  bimaterial.  By  using  an  entirely  different  methodologies, 
the  most  singular  solution  of  the  steady  state  problem  was  obtained  by  Yang  et  al. 
(1991)  and  the  complete  solution  of  the  steady  state  problem  was  given  by  DENG  (1992). 
However,  in  the  present  analysis  the  functions  Aj^(f)  and  B^(t)  are  allowed  to  be 
functions  of  time. 


3.2  Solutions  for  0m(*? 1,772,*)  and  *?2,  *)  for  m  =  2 


For  m  —  2,  the  equations  of  motion  (2.12)  are  coupled.  They  take  the  form, 

1  2u1/2  r  i 

■  . 

1  2t>1/2  f  1 

V’JUl  Oft,  0  +  ^^2,22  {r}uV2,t)  =  -^rIrTl{vl/2Go(^;0}(< 


(3.36) 


where  Fo{zi\  t )  and  Go{z»i  t)  correspond  to  the  solution  of  (2.12)  for  m  =  0  and  are  given 
by  equations  (3.32)  and  (3.33). 


In  order  to  obtain  the  next  most  singular  term  in  <h{V\iV2^)  and  tM’7i<r?2i*)»  we 
should  only  consider  the  most  singular  terms  in  Fo(zi\  t )  and  Go{z,',  t).  Therefore,  for  the 
material  above  the  interface, 


F0(zf,t)  =  a0{t)Ao(t)zf+'1  +  bo(t)Ao{t)zf  “  ' 
G0(z,;t)  =  Co(t)Ao(t)^il  +  <*0(t)A0(<)^_t< 


(3.37) 


where  A0(t)  =  A<,0)(0,  given  in  (3.35),  and 


a0(t)  - 


[(1  Fa])  -2qa.]  e‘ 

flD(v)  cosh  C7T 
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m 

Co  (0 

4o(t) 


[(1  +  <*\)  +  2T}Qa]e-tr 
fiD(v)  cosh  en 

[2 at  ~  f)(  1  -|-  Q,)]  e<y 

fiD(v)  cosh  or 

[2a;  + 17(1  +  0*)]  e~‘* 
piD( v)  cosh  or 


For  the  material  below  the  interface,  we  need  to  change  the  parameter  or  to  —or. 


Substituting  (3.37)  into  (3.36)  and  carrying  out  the  differentiation  with  respect  to 
time,  (3.36)  becomes 


and 


(*?i ,*?2,0  H — o^2’22  ~ 

af 

~of^Re  [^o(0ao(02r  ~  ^4o(<)6o(<)^r*e 

y/va, 


z?  In  zi 


2  at 


B0(t)a0(t)zr  +  B0(t)bo(t)zf 

■Jvai 


z(  3z, 


+  (\/uv4o(<)ao(<))  +  Bo(t)a0{t)  z{ 


+  IC 


+ 


\ji  (^WMO)  +  &>(«)M0  2*  j 


^2,11  +  *4  02  >22  (>?I,  >?2,0  = 

Qi 

2y/v 


ajcj 

y/vat 


^Im  [ikyfi  [+0(f)cb(<)<  -  A0(t)d0(t)zr]  zi  In ; 
-  [&(«)<*(*)**  +io(<)<*o(*)*rt]  z7*z, 

+  (v^o(0<*(0)  +  ~~£o(0«>(o]  +“ 


(3.38) 


> 


(3.39) 
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where 


^o(0  —  ( -  + ie j  A>(0  , 


Bo(t)  —  (  2  +  *e)  (2  +  zc )  ^o(0  • 


The  most  general  solutions  to  equations  (3.38)  and  (3.39)  are 

<h{ni,V2J)  =  Re{F2{z,;t)  -z,F{z,]t)  -  zjF(zt;t )} 
=  Im{62(2,;*)  -z,G(zt;t)  -l]G(z,\t)} 


where 


F(z-t)  =  D^aoitfiz^  +  DdboWz^-' 


and 


i-. 


>  • 


KUi) 


(3.40) 


+  e  {/v,(0oo(0^+,t  +  A'/(0M0^"‘}  In  ' 

G(z;t)  =  Da{co(t)}z}+'‘ +  D.{d0(t)}z}-“ 

+  e{Ka(t)co(t)z$+i‘  +  K,(t)d0(t)z}-'‘}lnz  J 

F(z;t)  =  ^(0ao(0^+U  +7?t(t)bo(t)si- 

G(z;t)  =  B,(t)co(t)z$+t' +  B,(t)d0(t)zl 
The  two  operators  /)/{■}  and  Dt{-}  are  given  by 

DMt)}  =  2^{^  +  'e)  ^  [V'/2p(<)  (I  +  k)  Mt). 

+  »2o>>  +  *e)  +  z  “,1/2p(0'4o(<)  j 

where  p(t)  is  a  real  function  of  time  t.  Also 

B‘M = (5 + i£) Mt)  | 

it?Ao(0 


In  (3.40),  F{zi\t),  F(zr,i),  G(z,;t),  and  G(z,;i)  are  totally  determined  by  F0{zf,t)  and 
G0(z,;t),  given  in  equation  (3.37).  The  coefficients  of  functions  F(z,;t),  F(zf,t ),  G(z,;t), 
and  G{z,\t )  are  related  to  the  crack-tip  acceleration,  the  time  derivative  of  Ao(t),  as  well 
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as  the  crack-tip  speed  and  Ao(f)  themselves  through  the  definitions  of  Di',{A0(t)} ,  Biit(t), 
and  Ki',(t).  It  should  be  noted  at  this  point  that  these  definitions  reduce  to  the  equivalent 
ones  corresponding  to  the  transient  crack  growth  in  homogeneous  materials.  Indeed,  if  < 
is  set  to  be  zero,  the  expressions  for  Di,{A0(t)}  and  Bit,(t)  that  appear  in  FREUND  and 
ROSAKIS  (1992)  are  obtained.  Once  again,  it  is  clear  that  for  the  steady  state  situation, 
functions  F(zy,t),  F(z/;t),  G(z,\t),  and  G(z,;t)  will  vanish.  The  undetermined  functions 
F2(zr,t)  and  G2(z,;t)  are  analytic  in  the  upper  half  plane  for  the  material  above  the 
interface,  and  in  the  lower  half  for  the  material  below  the  interface.  These  functions 
are  at  the  moment  unknown  and  will  be  determined  below  by  using  the  boundary  and 
bonding  conditions. 


Associated  with  <t>2{ViiV2^)  an^  the  components  of  displacement  will  be 

«!2)  =  Re  {/*(*;  O  +  a.^*.;*) 

-  [z/F'U,;  t)  +  zi]F'(zi;  t)  +  F(z, ;  <)  +  2 z,F{zi;  t)J 

-  as  [z,G'(z,;t)  4 -z2,G'{z,;t)  -  G(z,;t)  -  2z,G{z,;  t)j  J 

-  ,  (3-41) 

u[2)  =  Im  {a,F2'(z,;  t )  4-  G'2{zs;  t) 

-  at  [z/F'(zt;0  +  z^F'{zl;t)  -  F(zf,t )  -  2 z,F(z,;t)] 

-  [ zaG'(z, ;  t )  -f-  z2tG'(z,;  t )  4-  G(z,;  t)  4-  2z,G(z,;  <)]  | 

and  the  components  of  stress 

off  =  fiRe  {  (l  4-  2a?  -  a,2)  F?(zr,  t)  4-  2a.G"(z,;  t) 

-  (l  +  2af  -  a])  [z,F"(z,;  t )  +  z*F“(zr,  t)  +  2 F(z,;  <)] 

-  2  (l-^)  +  2Q'1-r-^2Q-  \F'(zr,t)A2zlF'(zrJ)\  ’  ^ 

-  2a,  \z,G"{z,\ t)  +  Ti]G"{z,-,t)  -  2 G(z,; t)] } 
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-  (l  +  a])  [*,F"(*i;  t)  +  z?F"(zr,  i)  +  2 F(z,;  <)] 


-  2 


O-^2) 


2  (af  —  a*) 


1  -or? 


[F'{z,]t)  +  2zlF'{zl;t)\ 


-  2a,  \z,G"{za- 1)  +  t)  -  2G(z ,■  <)] } 


(3.43) 


and 


(3.44) 


*i?  =  ~»lm  {2 oiFH(z,;  t)  +  (l  +  a])  G''{z$ ;  <) 

-  2a,  [z,F"(*,;  0  +  t)  -  2F(*,;  i)J 

“  (l  +  *2)  \z,G"(z,;t)  +  z]G"(z,-,t)  +  2G(z,;0] 

"  2  (l  -  a,)  [(?'(*,;  <)  +  2I,<7'(*,;  <)] } 

To  produce  a  more  compact  form  for  the  boundary  and  bonding  conditions,  one  needs  to 
define  the  following  quantities:  First  let  P *,  Qk ,  U k,  and  V\  be  obtained  from  matrices 
•Pfci  and  V*,  respectively,  by  changing  the  sign  of  the  off-diagonal  elements,  and 


let 


Mfc  = 


fim(v)  0 

Nk  = 

Hm(v) 

0 

0  /zn(n) 

k 

0 

-fin{v) 

where 

m(v)  =  (l  -  a?)  -  3 

Also,  define  complex  vectors, 


2  (qf  -  a?) 


n(n)  =  1  -  a]  . 


fk{z\t)  =  (  F2*(r;t),  G2*(z;f)  )T  > 
/*(*;»)  =  (  /ifel),  6,(r;i)  )T, 
/»(*;<)  =  (fi(r;<),  Gt(2;())T. 


Then,  by  using  the  above  definitions,  the  traction  free  condition  on  the  crack  faces  will 
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be 


Pi  [/i+(»?i;0  -  viK+  fat)  -  nl'fT  fat) 

~2 Mx  +  -2NX  fat)  +  2nxfx  fat) 


-2  Px  fifat)  ~2Qifi  fat)  =  o 


fT  fat)  -  m  fl  fa1)  -  tff"  0?r,*)] 

_ .  1*«+  _-TV/+  1 

fi  fat)  -mf 2  fat)-vlf2  (>h;*)] 

—2M 2 j/2  (»7i;*)  +  2»?1//3  (»?i; ^)]  -2iv2j/2  fat)  +  2mf2  fat)\ 


-2  p2  f2  fat)  -  2  Q2  f2  fa  t)  =  o 


The  continuity  of  traction  along  the  interface  will  reduce  to 

{■Pj  [f!+fat)  ~  Vi  fT  fat)  -  vlK+fat)] 

+Qi  ffi-fat) ~mh  fat)-vVfi  (»7i; <)] 

-2MX  [/i+(T?i5<)  +  27?1/'+(r;,;t)j  -2NX  fat)  +  2r)xfx  (t?x; «)] 
'2  Px  fifat)-2Qi  7T(*h;<)} 

-  {^2 \f2~fat)  -  mfi  fat)  -  »?i/T(»?i;*)] 

+Q2  [TT^iJO  -mf?  fat)  -Vih*  fa*)] 

~2M2  [ft  fat) +  2ih fifat)]  -2 N2  f2  fait)  +  2vi/2  (»7i; <)] 
-2  P2  f2fat)-2Q2  f2fat)}  =  o 


}iVm<  o. 


(3.45) 


,V  T}i  >  0, 


(3.46) 
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and  the  continuity  of  the  displacement  along  the  interface  will  be 

{Ur  -  vlK (*/»;<) 

+v:  [7'r(i?i; t)  -  -  v(f7 to; *)] 

-  Ui  +2vifUvi’t)}  -  Vi  [7i  fti;*)  +  2>?Ji  (71;*)]} 

-{e72  [/a'C’JliO  -  (wO-Vi/i  fo»;0] 

+v2  [?2+  (?ii;0  -  vif*  (*n;0  -  vff?  (w<)] 

-  LT2  +  V’a  [K(n»;0  +  2i?iK 


,Vr?,  >0. 


(3.47) 


Similar  to  the  procedure  in  the  previous  section,  by  rearranging  the  bonding  condition 
(3.46)  and  (3.47),  we  may  introduce  two  new  functions  n{z\t)  and  0(z-,t),  which  are 
analytic  in  the  cut-plane  S+  U  S~ .  In  order  to  keep  our  notation  short,  we  define  some 
new  quantities, 

gk(z;t)  =  -  zfl(z;t)  -  z2~f'l(z;t)  1  .  .  .. 

>  ,  k  €  iM). 

-  2P;lMk[f'k(z;t)  +  2zf'k(z;t)\  -  2P;1  Pk  ?*(z; 0  J 
Therefore,  we  can  write  that 

«(*;<)  *  Pi9i(*n)  -  1 


ff'(z;t)  =  t/i0j(z;t)  —  V?  ff2(z;f)  +  0i(*;O  — 


,  z  €  5+  ,  (3.48) 


k(z;<)  =  Pig2{z\t)  -  Qi0i(*;O 
0'(z;t)  =  UtfatoO  “ 


" 

■; t)-qx(z;t ) 


z  €  5-  ,  (3.49) 


where 


()  =  2  (I»M,  -  /)  [£(*;  0 + +  2  (**  -  n*)  f,(*i ») 

(*;<)  =  2  (i,  2V,  -  j)  [/!(*;()  +  2*X(*;«)]  +  2  (£»«,  -  V,)  ?»(*;«) 


(3.50) 


and 

I=\l  °1 
1  [  0  i  ’ 

Vector  qk  ( z;t )  is  related  to  vector  qk(z;t)  by 

kk  (z;t)  -  (  -qk]{z\t)  )T  ,  qk(z;t)  =  (  qk\z',t),  ql2)(z;t )  )T  . 

This  notation  will  be  used  throughout  the  paper  to  signify  this  operation.  In  calculating 
qk  ( z\t )  and  qk(z;t)  in  equation  (3.50),  we  have  used  the  fact  that 

PklMk  =  Q-k'Nk  ,  p;1  Pk=  Q r  Qk,  *  €  {1,2}. 


1  0 
0  -1 


By  solving  equations  (3.48)  and  (3.49)  for  gk{z',t)  and  gk{z]t),  we  obtain 
*(*;  t)  =  Px'H1  [e'(z;  t)-  L2  K(z;  t)  -  [9l(z;  t)  -  ~q2{z ;  <)]  } 

g2{z;  t)  =  Qi'H'1  { 6'{z;  t)  -  Lxk(z ;  *)  -  [fl,(z;  t)  -  q2(z;  <)j  j  J 


LzeS+,  (3.51) 


and 


g,frt)  - -P?  A'*  ■ 

6'{z-t)~  Lx  k (z;t)  - 

[«a(2;«)  —  9x(^;<)J 

}) 

j 

?,(*;<)  =  -Qr‘  Hl  ■ 

0'{z-t)  -  L2K(z;t )  - 

}J 

,2  €S~.  (3.52) 


It  can  be  seen  that  the  above  equations  are  very  similar  to  equations  (3.14)  and  (3.15) 
with  the  exception  of  terms  q\{z',t)  —  q2{z]t)  and  q2(z’,t )  —  <jh(z;t),  which  are  totally 
determined  by  the  solution  for  m  =  0. 


By  substituting  equations  (3.51)  and  (3.52)  into  the  boundary  condition  (3.45),  one 
can  show  that  k(z;  t)  is  an  entire  function.  As  a  result,  the  boundary  condition  (3.45) 
will  reduce  to 

==**(»?,;<) +  *(i?i;<) ,  v*,  <o  ,  (3.53) 

where 

#c(q,;f)  =H  [fli’faiiO-iafoi'.O]  “  [«2  ”  9i  fo;1)]  • 
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Equation  (3.53)  also  represents  a  Riemann- Hilbert  problem  for  0'(z;t).  It  requires 
that  0'(z\t)  is  analytic  in  the  cut-plane  S+  U  S~ ,  and  along  the  cut  satisfies  the  above 
equation.  By  using  the  properties  of  our  asymptotic  expansion,  equations  (2.9)  -  (2.11), 
it  can  be  shown  that  0’(z ;  t)  should  vary  as, 

|  |  =  O  (  \z\°  )  ,  as  \z\  -►  0  ,  (3.54) 


for  some  a  >  0.  The  complete  solution  of  (3.53)  is  generated  by  splitting  the  problem  to 
the  following  two  parts. 

°  > 

To  obtain  the  first  part,  let  6  {z-,t)  be  an  analytic  function  in  the  cut-plane  S+  U5~, 
such  that 


>'+ 


H0  (r]1-,t)-H8  fat)  =  RK{T)ut)  ,  V  T/i  <  0  . 


(3.55) 


This  is  exactly  the  same  as  equation  (3.18).  One  basic  difference,  however,  is  that  unlike 

o 1 

the  previous  case,  here  0  (z;  t)  has  to  satisfy  (3.54)  (recall  that  before,  a  >  -1).  As  a 
result  of  this  observation,  in  the  material  above  the  interface,  the  solution  for  gx  (  )  is 

given  by 


(1)r  [(1  +  Q* )  ~  2^o4)  etn  x+ft  o 

■  - JDM  cosh,,  *  A,{Z't] 

1(1  _+_gl.)  +  ,) 


fiD(v)  cosh  C7T 


and 


1  ll 

T  +  aJ 

2a,  \ 

»D{v)  \  \ 

kl  + 

1  +»,) 

(2) 

1  (*;<  = 

(2 Qi  -  r)(l  +  aj)]  e' 

..  T\t  _ t- _ 

(3.56) 


[2o,  +  »7(1  .  . 

+  iiD(v)coshtn  “ 


(3.57) 
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where 


o  .  .  (  °d)  .  .  0(2)  \T 

gl(z\t)=^gl  (z;t),  9 j  (*;<)]  , 

0  0 

and  the  entire  functions  Ai  {z\t)  and  B2  (z;<)  can  only  be  determined  by  the  far  field 
conditions.  Similarly,  the  solution  for  g2  (*;£),  in  the  material  below  the  interface,  can 
be  obtained  by  changing  the  corresponding  parameters  in  equations  (3.56)  and  (3.57). 

The  second  part  of  the  solution  is  obtained  by  letting 

0\z;t)  =  $'(z;t)~  0  [z\t) 

£k(z;*)  =  9k(z\t)~  9k  (*;<) 

Then,  6  (2;  t)  will  be  analytic  in  the  cut-plane  5+U5‘,  and  satisfy 

H  d'+{rji'>t)  -  HO'  (i?,;<)  =  ,  V  ^  <  0  .  (3.58) 

Because  the  right  hand  side  of  equation  (3.58),  k(rji;t)  is  totally  determined  by  the 
solutions  (f>o(r]i,r)2,t)  and  i/’ofaii 272,  t),  S'(z]t),  and  therefore,  gk(z;  t)  are  also  completely 
determined. 


By  using  the  results  in  Appendix  II,  we  can  write 


47Tt  /-00  I A2  L+{i}i)  M  L  (r)i)  )Vi-z 


where 


L(Z)  *  22+i 


The  explicit  dependence  of  qk(z-,t)  on  z  can  be  obtained  from  equation  (3.50), 
qk(z]t)  =  t*.{a0(<), co(<)}^*+ie 

+  c  (3  +  2it){LkMk  —  /)4*{«o(t), co(/)}2a+,t  In  2  ►  ,  (3.60) 

+  i(3-2i€)(LkMk-I)kk{bo(thMt)}z"-i<lnz  d 
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where  operators  like  tk{a0(t),  co(t)}  and  fc*{a0(t)>  <*>(<)}>  etc->  are  given  in  Appendix  I. 
From  the  definitions  for  qk  ( z ;  t)  and  the  above,  one  can  get 

ql(z;t)-^2(z;t)  =  /3z*+t'  +  i  £z±+it  In 2  -  i  h*~u  In 2  ' 

-  -  ,  (3-61) 

q2(z;  f)  —  q^iz;  t)  =  yz^+,t  —  f)z*~11  -f  c  C22+,t  In  2  —  i  £zi~u  In  2 
where  quantities  of  (3 ,  7,  and  «  are  also  given  in  Appendix  I.  It  should  be  noted  here 
that  /3  and  7  depend  on  the  crack-tip  speed  and  the  complex  parameter  Ao(t),  as  well 
as  their  time  derivatives.  However,  £  and  c  depend  only  on  the  crack-tip  speed  and  the 
complex  parameter  Ao(<)-  The  right  hand  side  of  equation  (3.58),  k(r)\\t),  becomes 


*0?i5<)  =  c  wd(-i7i)s+i<ln(-?71)-|-cu>j(-i7i)i~’tln(-71) 


,  V7l<0,  (3.62) 


where 


=  *  |e“£'  H  £  +  etirHs} 

*  . 

u;t  =  t  |e~4,r  H  f3  +  e<,rtf  7}  -  ire  { e’"  ff  (  -  e" JT«} 

Once  again,  it  can  be  seen  that  ujj  does  not  depend  on  the  time  derivatives  of  the 
complex  parameter  Ao(t)  and  the  crack-tip  speed,  while  u tt  depends  on  these  quantities. 
The  functions  inside  the  integrand  of  equation  (3.59)  can  be  rewritten  as 

.  u/  _  \  1  i  *.  <  „  2»t  i_/ _ \  ) 


~ =  -te,rt|cu;rfln(-i7i)-(-cu>d(-7j)~2,‘ln(-j/i) 

+u;t  +  w^-i?!)"2"} 

=  -te-ir<{c  ^(-^"lnf-ihj  +  e  <lrfln(~»7i) 
£  (»?i)  l 


_  • -»* 

^  Oft) 


(3.63) 


In  order  to  obtain  the  solution  for  ft (z\  t),  we  recast  equation  (3.59)  into  the  form  of  a 
Stieltjes  transform  by  using  (3.63).  However,  one  can  see  that  for  our  case,  a  closed  form 
evaluation  of  the  Stieltjes  transform  integral  is  very  difficult.  At  the  beginning  of  this 
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section,  it  has  been  mentioned  that  only  the  most  singular  terms  in  the  solution  of  F0(z;  t) 
and  Go(z;t)  are  considered.  This  implies  that  we  are  only  interested  in  the  region  where 
|z|  — ♦  0,  i.e.  very  close  to  the  interfacial  crack-tip.  As  a  result,  instead  of  evaluating 
the  entire  Stieltjes  transform,  we  only  need  to  study  the  asymptotic  behavior  of  that 
transform  as  \z\  — ►  0.  The  details  of  this  asymptotic  analysis  are  given  in  Appendix 
III.  If  only  the  leading  terms  in  the  Stieltjes  transform  are  retained,  by  using  the  results 

a/ 

provided  in  Appendix  III,  the  solution  for  6  (z;t)  can  be  obtained  as 

0(z;t)  = 

+ 

where  in  developing  the  above  equation,  the  relation 


- 


(In  z)^  + 


-  £ tzs 


In  z 


Cttz2+U  ~  Cuz: 


+  0(1*1) 


(3.64) 


A!  ’ 


has  been  used,  and  the  following  notations  have  been  defined: 


Cd  "  ^T,rUd 


C.  = 


e-**  f 
4ttA,  \ 


r  u>t  — 7 


tnt 


sinh(27rc) 

<■  ■  £{[?'— 


|Vecosh(27re)  .  iw 

+  — rrrrz — r~  •*  -t-  —  .  — — 


sinh2(27rc) 


sinh(27rc) 


HI  j 


In  constructing  the  entire  solution  for  0,(z;t)  and  g2(z;t),  the  leading  terms  in  (3.56) 
and  (3.57),  are  considered.  This  is  consistent  with  the  fact  that  (3.64)  contains  only 


33 


leading  terms  of  the  same  order.  The  final  solutions  for  gx(z\  t)  and  g7(z]t),  are  therefore 


-  qt'  h'  cj2i-" 


(lnz)2 


+ 


+ 


P;'H-'  «,  -  i  ()  zi+“  -  Q-‘  h'  (<,  -  i  t)  zi- 

«„  h'1  (c„  -7) 


In  z 


+  S*i+“  A  w  +  + 1 0  ( w  > 


J ,  z  6  5+,  (3.65) 


and 


9i(*;  0  -  -< 


Pj‘  H  0*’s+"  -  Q;'H-'(dzt 


.  -1 


4+«£  _  0-1  W-V  ,o-«« 


(lnz)2 


P?1  ff"‘  «,  -  i  t)  ^+"  -  (f,  -  i  ()  »!-” 

PJ1  h”  ((„  -  -r)  2i+”  -  Qi'H~'  ({„  -  3) 


In  z 


«"Qr<  i-.T 


y,2£5-,  (3.66) 


+ -^hirz?+“  a  w + + o  <  w ) 

where  A2  (t)  =  A2(<M)- 


Our  final  target  is  to  find  the  complex  potentials  /*(z;f)  =  (  F2k(z’,t),G2k{z\t)  )T  ■ 
After  some  manipulations,  f"(z;t)  and  /^(z;  t)  can  be  expressed  as 

/?(*;  0  =  i  |prlH-‘<^>+,<  -  Q;'  h'  (In  zf 

+  |  [PJ-'ff-1  (C.  -  i  0  +  i  «-«{oo(i),  co(t)}]  ^+“ 

+  |  [P;'H~'  «„  -  0)  +  »nW(),  <*>(!)}]  I’+,< 

-  [<?:’  »  '  (?k-7)-®<i{M0.<«01  *H'} 

^-"Aa(()+0(W) 


etnP71C  1  o  1  ” 

+  -  V  C*>+“  /t,  (<)  + - 

cosh  tir  cosh  cir 


,  z  €  S+,  (3.67) 
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and 


(In  2)' 


*+•< 


In  2 


1 ,  2  €  S~,  (3.68) 


zi+<’ 


fUz;t)  =  — «  |fj‘  h'‘  C„zJ+"  -  Q?H-‘(dzi-’ | 

((,  -  e  «)  -  e  u>,a(iii)(<),co(()} 

-  [Qj-'H"1  (C,  -  i  I)  +  i  ®«{M0.  *(<)}]  **■"] 

(C«-7)-®aW(),co(()) 

-  [Qj-'h-1  (C„  -3)  +®o{W0. *(«))] 

+  A,  (()  +  ^Xh<A,(()  +  O  (  W  ) 

where  the  operators  «;,**(•,  •)  and  «>**(•,•)  are  given  in  Appendix  I.  By  integrating  the 
above  two  expressions  with  respect  to  the  complex  argument  z ,  we  can  finally  obtain  the 
complex  potential  fk(z ;  t )  =  (  F2k(z;  t),  G2k(z‘,  0  )T  for  both  materials. 


Since  equations  (3.67)  and  (3.68)  are  directly  related  to  the  stress  components  around 
the  interfacial  crack-tip,  some  of  the  noteworthy  features  of  the  asymptotic  field  can  be 
studied  through  them.  The  most  interesting  feature  is  that  there  exist  two  terms  in  the 
above  equations,  which  are  totally  different  in  nature  from  the  terms  found  in  the  solution 
of  a  crack  propagating  transiently  in  a  homogeneous  material.  The  first  of  these  terms 
is  that  associated  with  z^^ln  z )2.  This  term  is  clearly  associated  with  the  interfacial 
nature  of  crack  growth  since  it  is  proportional  to  the  quantity  c.  This  quantity  is  also 
related  to  the  transient  nature  (existence  of  non-zero  accelerations)  of  the  problem.  By 
observing  that 

dc  dt> 
dv  dt 

One  can  easily  see  that  i  and  thus  the  z1/r2(lnz)2  term  vanish  either  when  the  crack-tip 
speed  is  constant  and/or  when  the  material  mismatch  parameter  c  vanishes. 

The  second  term  is  that  associated  with  z1/2  In  z.  The  coefficient  of  this  term  is  related 
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to  the  complex  parameter  Ao(t)  and  also  depends  on  the  crack-tip  speed,  as  well  as  on 
their  time  derivatives.  So  it  depends  on  both  0(1 )  and  Kd(t).  It  can  be  seen  that  for 
constant  speed  transient  crack  growth  ( v  —  0,  Kj  ^  0),  this  term  will  still  be  present. 
Indeed  the  r!^2lnr  term  has  been  observed  by  WILLIS  (1973),  who  studied  a  particular 
constant  velocity,  transient  interfacial  crack  growth  problem.  Both  of  these  two  terms 
which  include  logarithms  will  vanish  at  the  same  time  only  if  the  situation  is  strictly 
steady  state.  Otherwise  one  or  both  will  be  present.  These  logarithmic  sigularities  are 
the  consequences  of  the  existence  of  both  the  interface  and  the  transient  nature  of  the 
propagating  crack.  For  the  case  of  crack  growth  in  a  homogeneous  material  (c  =  0). 
(3  =  7  and  £  =  see  Appendix  I.  This  is  true  even  if  crack  propagation  is  transient.  As 
a  result,  it  can  be  shown  that  u><i  =  o>t  =  o,  and  consequently,  £t,  and  £t(  will  vanish. 
The  logarithmic  terms  also  disappear.  In  this  case,  the  transient  field  reduces  to  the  one 
obtained  by  LlU  and  ROSAKIS  (1992)  which  does  not  feature  any  logarithms.  It  should  be 
stated  at  this  point  that  transient  higher  order  terms  involving  logarithmic  singularities 
have  also  been  observed  in  the  solution  of  dislocation  lines  propagating  transiently  in 
elastic  solids  (Callias  et  al .,  1990,  and  MaRKENSCOFF  and  Nl,  1990).  These  terms 
were  shown  to  vanish  when  the  dislocations  propagated  with  constant  speed. 

In  this  section,  we  have  provided  a  procedure  which  allows  us  to  investigate  higher 
order  transient  effects  systematically.  By  imposing  the  boundary  and  bonding  conditions 
on  the  complex  potentials,  the  problem  was  recasted  into  the  Riemann-Hilbert  problem. 
By  solving  the  Riemann-Hilbert  equations,  and  by  evaluating  the  Stieltjes  transforms, 
the  higher  order  terms  were  obtained.  This  procedure  can  be  repeated  to  any  order,  and 
we  may  therefore  claim  that  it  is  constructive. 
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4  The  asymptotic  elastodynamic  field  around  a 
non-uniformly  propagating  interfacial 
crack- tip 

For  planar  deformation  of  a  homogeneous,  isotropic,  linearly  elastic  material,  the  ordered 
array  [ua>£Qp,<Ta0),  a,  /?  €  {1,2],  is  said  to  be  an  elastodynamic  state  in  the  absence  of 
body  force  density,  if  the  following  conditions  are  satisfied 

_  1  ,  ,1 

£<30  —  2  V^Ori 0  ~^^01O  ) 

<ra0  =  2p£Q0  +  \£yySa0  ’  a,/?  6  {1,2},  (4.1) 

ara0'0  —  pUa  J 

where  p  is  the  mass  density  and  A,  p  are  Lame  constants  of  the  material.  In  addition, 
the  field  quantities  uQ,  eQ<3,  and  aa0  must  satisfy  the  smoothness  requirements  outlined 
in  Wheeler  and  Sternberg  (1968). 

In  the  Cartesian  coordinate  system  (£1,62),  let  i,&2,0  and  0m(fi,£2,O  be  solu¬ 
tions  of  equations  (2.12),  m  =  0, 1, 2,  •  •  •,  such  that 

,  £2,  t )  _ ^  „ 

0m(6,  &,t) 

V,m4-n(^l ,  ^2,  t )  q 

for  any  positive  integer  n.  Thus,  1,^2, 0  and  0m(6>£ 2,0  will  be  two  asymptotic 
sequences  as  r  =  (£j*  +  £f)l/l2  — * ►  0.  Define  <^(6,^2, 0  and  0(£i,&,f)  by 

=  £*.«■  1.6.0 

m=0 

m=0 

Then,  the  array  \uo,£a0,(Ta0},  a,  ft  6  {1,2},  will  constitute  an  asxjmptotic  elastodynamic 
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state  as  r  =  (^  +  £\)1^  — ►  0,  if  it  satisfies 

=  4>,a  +ea04>,3 

£a0  =  |  (Ua,0  +  U/3,a  )  ►  ,  0,0  €{1,2}.  (4.4) 

<?a/3  —  2 [A£&0  +  j 

Now,  consider  a  planar  body  composed  of  two  homogeneous,  isotropic,  and  linearly 
elastic  materials  bonded  along  a  straight  interface.  Let  the  two  displacement  potentials 
for  each  material  be  given  by  (4.3),  where  each  term  of  the  asymptotic  series  is  the 
solution  which  has  been  discussed  in  the  previous  section.  The  asymptotic  elastodynamic 
state  near  the  non-uniformly  propagating  interfacial  crack-tip  can  then  be  obtained  from 
relations  (4.4). 

For  its  importance  in  the  experimental  investigation  described  in  section  7,  we  only 
provide  the  asymptotic  expression  of  the  first  stress  invariant  around  the  interfacial  crack- 
tip.  However,  in  order  to  shorten  our  expression,  some  notation  needs  to  be  defined  first. 
In  the  expressions  below,  the  superscript  (1)  or  (2)  denotes  the  components  of  the  vectors 
defined  in  Appendix  I  and  in  previous  sections.  For  the  material  above  the  interface,  we 
may  define  the  following  quantities, 

n i(t)  =  -  '  iA;  {  |(1  +  a))hn  +  2o.ft,i)  -  [(1  +  ol)hl2  +  2  a.*,,)  iff1}  , 

Od  (*)  =  A,  {  [(1  +  Q'^u  j  d *  +  [U  +  +  2Q.hu]  cl'}  . 

( t(1 + a-)h" + 2a-*>'l  (<*"  - 4  fl”) 

-  [(1  +  a])kn  +  2q,Au]  (Ct(2)  -  e  f(2))}  4-  c  «□(*)}  - 

«■<*>  =  + +  2o-,‘2']  (^”  - i 

+  [{1+  a ])h12  +  2<a„]  (?!2)  -  e  ?<2))}  -  <  - 
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—  [(1  4-  <*l)hi2  +  2a3hn]  (Ct^  ~  /^2*)j  +  co(0}  > 

n,,(i>  =  ( t(I  +  Q-)h"  +  2a,Al’l  (^"  -  7<”) 

+  [(1  4-  ot2,)hi2  4-  2Q4/inj  (cL1  —  7*2)^  |  —  o(t)}  • 


Now,  one  can  show  that  by  using  equations  (3.32)  and  (3.67),  the  first  stress  invariant  in 
the  material  above  the  interface  will  be  given  by 


a\\  +  <722 

2/i  (a?  -  a?) 


Re  |  Ao  ( t)a0{t)zl  2+“  4-  Ao{t)bo{t)zJ J~‘e 


^ fiD(v)(l+  u  )  041  ^  +  Al^')  +  (*)a°(*)*2  +  + 


+ 


Ax  4-  Bit(t)ziz^11  +  Cu(t)zfz^+,t 


-  An  ( t)zf~'1-  Bu  (t)z,z;l'~u-  Cu  (t)zu;^] 

+  j4t(t)2/J+,e  4-  Bt(t)zizJ~3+,t—  At  (t)zf-'1-  Bt  (0*, ln2:' 
4-  i  [f h(t)z?+u-  Ud  (<)*/"']  (ln^)2}  +  O  (  \zi\  ) 


} 


(4.5) 


where 

A„{t)  =  n„(t)  -  {(|  +  it)  A!-o(‘)}  +  i  WWO}  -  2B, (()<.„(!)  , 

A„  (i)  =  h„  (()  +  ~  «) *.<««)>  +  f 7?, (<)«<)}  +  2S.COM.)  , 

B«(()  =  -  ( |  +  it)  (j  +  it)  AMO)  -  2f(l  +  it)K((i)oo(l) 

Bn  (<)  -  -  it)  (|  -  it)  ACWO)  +  2f(l  -  ie)7f,(0W0 
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(3  - !£)  BMb°w  ■ 

C„m  =  (i  +  tJ)  B,(t)a0(t)  ,  c,,(t)  =  -(j  +  <!)  W)WO  , 

A,«)  =  «■(<)- ^rfp  (|  +  «)  A-,(i)a.W  , 

A. w  =  n,  («)  +  2(‘  +  °|)c  (|  -  .<)  T,(«)M<) , 

B. (<)  =  -i  (|  +  «)  (5  +  «)  W)°o (0  , 

B.{t)  =  i(|-ie)(5-<t)^<(0M0, 

_ .  » _ m 

"  (f^KH) ' 


00  o 

In  expression  (4.5),  functions  of  time  Ao  (*),  Ai  (<)>  and  ->42  (<)  are  undetermined  by 
the  asymptotic  analysis.  On  the  other  hand,  functions  -4tt(t),  Cu{t),  Att  (<)>  ‘ '  •> 

O 

are  known  in  terms  of  Ao  (*),  the  crack-tip  acceleration  and  the  time  derivative 

O 

of  Ao  (t).  As  a  result,  these  functions  are  also  undetermined  by  the  asymptotic  anal- 

O 

ysis.  However,  their  dependence  on  time  derivatives  of  v(t)  and  Ao  (f)  constitutes  the 
mathematical  demonstration  of  transient  effects. 

It  is  often  convenient  to  express  the  first  stress  invariant  in  terms  of  real  quantities. 
For  any  complex  function  of  time  W(t),  let  its  magnitude  be  denoted  by  |VFj,  and  its 
phase  be  denoted  by  $(  W).  Meanwhile,  a  scaled  polar  coordinate  system  (rt,  $1)  centered 
at  the  moving  crack-tip  is  defined  by 

H  =  {tf  +  • 

The  first  stress  invariant  in  the  material  above  the  interface  can  therefore  be  expressed 


40 


as 


<?n  +  °22 


2 fi  (af-^of)  ~  I  Ao  (01  {So(^)c°s(clnr/)+  E0  (0/)  sin(e  In  r,)j  rf 

4a, 


1/2 


+ 


nD(v){\  +  a>4) 


|  A i  (<)|cos$(,4i) 


4-  t  |Sd(0/)cos(einr/)4-  Ed  (6t)  sin(c  In  r,)|  r//2(lnr/)2 
+  jEj(0/)cos(elnrj)4-  Si  (0/)  sin(elnr/)|  r}/2  \nr, 

4-  | Ett(0/)  cos(e In  r/)+  E«  (0/)  sin(e  in  r,)j  r//2 
+  [  A2  (01  |s2(0i)cos(clnr,)+  £2  (0/)sin(eln  r,)j  r]/2 

+  0(r,) 


(4.6) 


where 

Eo(«i)  = 

Eo  (0i)  = 

S2(0/)  = 
E2(0/)  = 
£<*(*,)  = 
E</  (0/ )  = 

St(0i)  = 

+ 


)■ 

)• 

-  |a0(Oe_<^  sin  4-  $(A2)j  -  MO**1  sin  ^  -  $(A0)j  J  , 

cos  ^  - 1  n d  (Ole'S‘ cos  +  $(fldYj  , 

-  ||nrf(0|c-ttf,sin  ^  +  |  hd  (01^  sin  ^  4-  $(S7d)^  j  , 

\M*) \€~t$‘ cos  -  I  At  (Ole'*'  COS  ^  4-  $(At)j 

|Bt(0|e“'*'cos  -  $(St)j  -  |  Bt  (<)|e'®'cos  -  $(B,)^ 

+  $(ftd)^  | 


2e  ||nrf(<)|e  t6‘  sin  j 'j  +  $(fid)^  -  |  hd  (<)(ce®*  sin 


ao(0e~tfll  cos  ^  -  $(Ao)^  +  bo(t)et6‘  cos  +  $(Ao) 

/ e .  o  ' 

^  +  *(  A>) 

a0(t)e~tdl  cos  ^  +  $(Aa)  j  +  bo(t)etB‘  cos  ^  -  $(A2) 


a0{t)e 


sm 


(f  -  *(A„>)  - 


bo(t)e 


tBi 


sin 
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E*(ft)  = 


+ 

s«(ft)  = 
+ 

+ 


+ 


ZttW  = 


+ 

+ 


+ 


-  sin  +  *(  A,) )  +  |  A,  <i)|e‘#'  sin  +  *(  A,)U 

sin  (y  -  +  |  Bt  (Ole"'  sin  ^  -  4>(B,)) 

(f 

cos  ^  -  |  Att  (0IC<#‘  cos  +  $(A(()j 

\Bu{t)\e~t9‘  cos  -  |  Btt  (Olc<e'cos 

|C«(<)|e_£tf'  cos  -  $(Cti)^  -  |  Ctt  (0le'*' cos  (^y  ~  *(C«)1 
sin  +  *(At)j  -  |  At  (*)|e*  sin  +  $( A, A  1  ft 

sin  -  $(£<))  -  |  Bt  (t)\et6t  sin  ^  -  *(£,))  J  0, 

i  ||firf(<)|C^  cos  +  $(0,))  -  I  nd  (t) |e<fi‘  cos  +  *(&))  |  ft2, 

-  1 1 Att(t)\e~*  sin  +  |  Au  (t)\e'*‘  sin  +  *(  A«))  } 

\Bu(t)\e~t91  sin  -  $(Bti)J  +  |  Btt  (0 \et6'  sin 

\Cu(t)\e-+  sin  -  *(Clt))  +  |  Ctt  (Ok*  sin  ^  -  *«?„)) 

{l^Wle^  cos  4-  *M,))  +  |  At  (t)\e<*>  cos  +  $(At))  }  ft 

cos  -  $(Bt)  j  +  |  Bt  (t)|e**‘  cos  ^  -  ${Bt)^J  j  ft 

c  sin  (^~  +  $(Hi) j  +  |  f U  (ON**'  sin  ^  +  $(Oi)^  |  ft2. 


+  $(Qd)  >  ft 


2c  lin^Ole-^  cos  +  $(fi^  +  |  hd  (t)\e‘*‘  cos 


The  first  term  in  equation  (4.6)  has  a  square  root  singuarity  and  oscillatory  nature. 


It  is  associated  with  the  complex  dynamic  stress  intensity  factor  Kd{t)  (defined  by  Yang 

O 

et  ai,  1991)  which  is  related  to  the  complex  coefficient  Ao  (t)  by 

A'd(t)  =  -2v^Ao  (0- 

The  second  term  is  the  so-called  T- stress  term,  and  is  independent  of  position.  The  first 
two  terms  have  the  same  spatial  form  as  those  obtained  under  steady  state  conditions  by 
DENG  (1992).  However,  the  remaining  four  terms,  proportional  to  the  square  root  of  the 
radial  distance  from  the  crack-tip,  are  more  complicated  and  have  some  unusual  features. 

O 

The  part  associated  with  j  Ai  ^)|  has  the  same  form  as  that  predicted  by  the  steady  state 
solution  and  is  of  order  r1/2.  The  term  of  order  r1/,2(lnr)2  has  a  coefficient  proportional 
to  e  —  e'(v)v(t).  This  term  vanishes  either  when  v  =  0  and/or  e  =  0.  The  remaining  two 
terms  contain  the  functions  £t(0/),  Sj  ($i),  £«(0j),  and  Ett  (0t)  which  depend  on  the  time 
derivatives  of  the  complex  dynamic  stress  intensity  factor  and  the  crack-tip  speed,  i.e. 
they  depend  on  transient  effects.  These  parts  also  vanish  for  steady  state  crack  growth. 
The  term  of  order  r1^2  lnr  was  first  observed  by  WlLLIS  (1973)  who  analyzed  the  stresses 
in  the  case  of  constant  speed,  transient  interfacial  crack  growth.  In  this  case,  v  =  0, 
Kd  0,  and  the  only  surviving  terms  will  be  of  order  r1^2  In  r  and  r1^2.  If  the  two  elastic 
materials  that  constitute  the  bimaterial  system  become  identical,  the  terms  associated 
with  r1/,2(lnr)2  and  r^2lnr  will  disappear.  However,  in  this  case,  the  functions  £«(0/) 
and  ( 6i )  do  not  vanish  and  reduce  to  the  ordinary  transient  term  given  by  LlU  and 
ROSAKIS  (1992)  in  studying  the  transient  growth  of  a  crack  in  homogeneous  materials. 
It  is  significant  to  note  at  this  point  that  transient  effects  may  noticeably  change  the  r 
and  9  structure  of  the  field  from  that  predicted  by  the  steady  state  approximation  (e.g. 
existence  of  logarithmic  r^2lnr  and  r1/l2(lnr)2  terms). 
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5  Properties  of  the  mismatch  parameters  in  dy¬ 
namic  inter  facial  fracture 


In  the  analysis  of  an  interfacial  crack  dynamically  propagating  along  the  interface,  there 
are  two  mismatch  parameters  which  depend  not  only  on  the  properties  of  the  materials 
that  constitute  the  bimaterial  system,  but  also  on  the  crack-tip  velocity.  The  properties 
of  these  parameters  are  very  important  since  we  have  seen  that  the  asymptotic  represen¬ 
tation  of  the  crack-tip  field  is  drastically  changed  due  to  their  presence.  One  of  these 
parameters  is  defined  by 


_  J_.  1^2  a  hu 

e  2t  n  1  +  p  ’  p  vtt  ’ 


(5-1) 


while  the  other  one  by 


In  the  above  two  definitions, 


(5.2) 


where 


hn  = 
hu  = 

h2i  = 


[2q7 a,  -  (1  +a?)l  f  2 aiat  -  (1  -1-  a])] 

l  t*D{v)  J ,  \  fiD(v)  J2 

fat(l  -a?)l  [0,(1  -a?)l 

1  M&(v)  ii  \  M&(v)  1 2 

,  f 01(1-02)1 

l  VD(V)  Ji  l  pD(v)  )2 


(5.3) 


(2  \  1/2  /  2\  1/2 

1  “  “2  J  »  Q*  =  V  Cj  /  ’  =  4a/°*  ~  ^  +  °*) 


To  illustrate  the  properties  of  the  mismatch  parameters,  we  choose  a  bimaterial  system 
composed  of  PMMA  and  AISI  4340  steel.  We  denote  PMMA  as  material-1,  and  AISI 
4340  steel  as  material-2.  The  mechanical  properties  for  these  two  materials  are  listed  in 
the  table  below. 
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Table:  Properties  of  selected  materials1 


Parameter 

fi(GPa) 

V 

ci(m/s)“ 

ci(m/ s)" 

c,(m/s) 

cR{m/s) 

p{kg/m 3) 

PMMA 

1.20 

0.35 

2081.7 

1761.5 

1004.0 

937.8 

1190.0 

AISI  4340 

80.0 

0.30 

5978.8 

5401.9 

3195.8 

2959.8 

7833.0 

'plane-strain  ''plane-stress 


For  both  plane  strain  and  plane  stress,  Figure  2  presents  the  variation  of  the  parameter 
T)  with  respect  to  the  crack-tip  speed.  We  can  see  that  r]  varies  smoothly  from  1.0  for  the 
stationary  interfacial  crack,  to  oo  as  the  crack-tip  speed  approaches  the  shear  wave  speed 
of  PMMA  (cW).  However,  the  situation  is  different  for  the  parameter  0.  In  Figure  3,  we 
can  observe  that  if  the  crack- tip  velocity  is  less  than  the  Rayleigh  wave  speed  of  PMMA 
(cg>),  0  varies  smoothly  and  tends  to  —1  when  the  crack-tip  speed  is  very  close  to  Cr*. 
Since  2?i(v)  will  change  sign  as  the  crack-tip  speed  crosses  Cr\  0  jumps  from  —1  to  1, 
and  then  tends  to  oo  as  the  crack-tip  speed  tends  to  c[lK  Figure  4  shows  the  variation 
of  0  when  the  crack-tip  speed  is  bigger  than  Cr\  Figure  5  presents  the  behavioT  of  the 
parameter  e  when  the  crack-tip  speed  is  below  the  Rayleigh  wave  speed  of  PMMA.  It 
shows  that  e  tends  to  oo  as  the  crack-tip  speed  is  very  close  to  Cr\  However,  as  the 
crack-tip  speed  crosses  the  speed  Cr\  since  0  is  larger  than  1,  c  will  become  complex, 
and  thus  e  can  be  written  as 

.  1.  -  1.0-1  .s 

£=f+2!'  '=2?b^TT-  (5'4) 


Figure  6  gives  the  variation  of  the  real  part  of  e  (i.e.  c)  with  repect  to  the  crack-tip 
speed  when  the  interfacial  crack  is  running  at  speeds  between  Cr*  and  c*1*.  We  can  see 
that  the  real  part  of  t  changes  from  —  oo  to  0  when  the  crack-tip  speed  is  in  the  range  of 


4J>  <  V  <  c<‘>. 


lthe  parameters  for  PMMA  are  from  CYRO  Industries,  Woodcliff  Lake,  NJ  06675;  the  parameters 
for  AISI  4340  steel  are  from  Aerospace  Structural  Metals  Handbook,  Battelle  Columbus  Laboratories, 
Columbus,  Ohio 
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6  The  asymptotic  field  of  an  interfacial  crack  prop¬ 
agating  at  a  speed  between  the  lower  Rayleigh 
and  shear  wave  speeds 


In  recent  experimental  investigations,  described  in  section  7,  bimaterial  specimens  com¬ 
posed  of  PMMA  and  AISI  4340  steel  have  been  tested  dynamically.  This  bimaterial 
combination  exhibits  a  remarkable  stiffness  mismatch.  It  was  observed  that  under  im¬ 
pact  loading  conditions,  interfacial  cracks  may  propagate  at  speeds  exceeding  see 
section  7.  This  experimental  observation  motivates  our  attempt  to  investigate  dynamic 
crack  growth  in  interfaces  at  speeds  exceeding  the  lower  Rayleigh  wave  speed.  In  homo¬ 
geneous  materials,  an  infinite  amount  of  energy  has  to  be  transmitted  to  the  crack-tip  to 
maintain  extension  at  the  Rayleigh  wave  speed  if  the  dynamic  stress  intensity  factor  is 
non-zero  (FREUND,  1990).  This  makes  it  impossible  for  a  crack  in  a  homogeneous  solid 
to  exceed  the  Rayleigh  wave  speed  of  that  material.  However,  for  a  crack  growing  along  a 
bimaterial  interface,  it  has  been  shown  that  as  the  crack-tip  speed  approaches  the  lower 
Rayleigh  wave  speed,  say  c^,  only  a  finite  amount  of  energy  has  to  be  transmitted  to 
the  crack-tip  if  the  dynamic  stress  intensity  factor  is  non-zero  (see  Yang  et  a/.,  1991). 
Accordingly,  there  is  no  energetic  restriction  for  an  interfacial  crack  to  exceed  the  lower 
Rayleigh  wave  speed.  Indeed,  the  experimentally  obtained  velocity  histories  reported  in 
section  7,  see  Figure  14,  are  seen  to  largely  exceed  the  Rayleigh  wave  speed  of  PMMA. 

In  the  analysis  of  previous  sections,  the  governing  equations  hold  for  crack-tip  speeds 
in  the  range  0  <  v  <  c^,  if  material-1  is  more  compliant  than  material-2.  Also,  the 
development  of  the  asymptotic  stress  field  around  the  tip  of  a  non-uniformly  propagating 
interfacial  crack  is  dependent  on  the  complete  solution  of  the  Riemann-Hilbert  problem. 
However,  from  the  procedure  provided  in  Appendix  II,  we  can  see  that  there  are  no 
restrictions  imposed  on  crack- tip  speed  from  this  procedure.  The  only  consequence  of  the 
restriction  that  the  crack-tip  speed  is  in  the  range  of  0  <  v  <  c^1},  is  that  all  parameters 
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appearing  in  the  solution  are  real.  Nevertheless,  the  mathematical  approach  is  not  limited 
by  this  restriction,  even  if  some  of  the  parameters  become  complex.  Therefore,  we  can 
directly  extend  our  solution  to  the  case  where  the  crack-tip  speed  exceeds  the  lower 
Rayleigh  wave  speed. 

Suppose  the  properties  of  the  materials  constituting  the  interface  are  such  that  c*,}  < 
cg\  and  c^  <  v  <  c^.  As  we  have  shown  in  the  previous  section,  the  parameter 
f?  remains  real,  but  e  becomes  complex  and  is  given  by  equation  (5.4).  If  only  the 
leading  term  is  considered,  under  the  requirement  of  bounded  displacement,  or  integrable 
mechanical  energy  density  (FREUND,  1990),  the  two  complex  displacement  potentials  in 
equation  (3.32)  for  the  material  above  the  interface,  become 


F0(zf,t) 


[(1  +  or*)  -2Tjaa)et 


(2  +  i  e)(l  +  t  e)nD(v)  sinh  e  x 
[(1  +<**)  + 2i?a#]  e~ 


,.2-f  te 


Mt) 


Go{za;t )  — 


(2  -  i  €)(1  -  i  e)(iD(v)  sinh  e  x 
[2a,  -  V(l  +  a2t)]e‘* 


(2  +  i  c)(l  +  i  e)fiD(v)  sinh  c  x 
[2a, +  *?(!+ a?)] 


—z2rlMt) 


(2  —  i  e)(l  —  i  e)/*Z?(t;)sinh  c  x 


(6.1) 


for  an  arbitrary  complex  function  Ao(t).  To  obtain  this  result,  the  definition  of  t  in  the 
speed  range  <  v  <  equation  (5.4),  has  been  used.  For  the  material  below  the 
interface,  we  need  to  change  e  x  to  -  e  x  in  the  above  expressions.  By  setting 


A0(t)  =  A(<)e'*<‘>  , 


the  first  invariant  of  stress  for  the  material  above  the  interface  becomes 
<7\l 


-f  022  =  - a> -  ^  ~  {27a,  cosh  [e  (x  -  0/)] 

D(u)sinh  t  it 


-  ( 1  -f  or?)  sinh  [e  (it  -  0/)] }  cos  In  r,  +  $(<)) 


1 


(6.2) 
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It  can  be  observed  that  oscillations  still  exist  along  the  radial  direction.  However,  there 
is  no  singularity  at  the  propagating  crack-tip. 


At  a  position,  r,  ahead  of  the  interfacial  crack-tip,  the  traction  on  the  interface  can 
be  expressed  as 

»?cr22(r;  t)  +  i<rl2{r;  t)  =  A0(t)  .  (6.3) 

At  a  position,  r,  behind  the  interfacial  crack-tip,  the  crack  face  displacement  difference 
is  found  to  be 


C  ,  Mr;t)  2 r,hl2  r1+,£  ,  ,.s 
«i(r;0  - »-  —  =  — r — -M*)  ■ 
v  sinh  e  1  +  i  t 


(6.4) 


If  the  interfacial  crack  extended  an  amount  6,  then  the  energy  released  by  this  extension, 
AH7 (5)  can  be  calculated  by 

&W(6)  =  lf‘  ;()<,(«  -  f,;0  +  <r,j({i;  <)«,(«  -  £ii<)}  •  (6.5) 


By  using  (6.3)  and  (6.4),  we  can  express  the  above  equation  as 

2rf2h12\Ac 


AW($)  = 


sinh 


e  [Jo  1  +  i  t  J 


(6.6) 


Further,  it  can  be  shown  that 

r‘  (i-J.),+i,fr” 


6  If,  _  c 


JO  1  -f-  t  £  ■'0  J  —  i  £ 

Therefore,  the  energy  release  rate  at  the  tip  of  an  interfacial  crack  moving  at  speeds  in 
the  range  <  v  <  Q,  will  be 


o  = 

6—0  6 


(6.7) 


This  result  may  be  anticipated  since  in  this  range  of  speeds,  both  stress  and  strain  are 
bounded.  Equation  (6.7)  states  that  if  the  speed  of  the  interfacial  crack  is  in  the  range 
c^  <  v  <  no  energy  is  needed  to  create  new  surfaces. 
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7  Experimental  evidence  for  the  importance  of 
transient  effects  in  the  dynamic  fracture  of 
bimaterials 


To  investigate  the  validity  of  the  analysis  presented  in  this  work,  a  sequence  of  dynamic 
impact  experiments  of  bimaterial  specimens  has  been  performed.  Stress  waves  generated 
by  impact,  load  an  interfacial  pre-crack,  which  subsequently  propagates  dynamically 
along  the  bimaterial  interface.  High  speed  interferograms  of  the  near-tip  region  of  the 
propagating  crack  are  recorded.  The  optical  method  used  is  the  newly  developed  method 
of  Coherent  Gradient  Sensing  (CGS)  (TlPPUR  et  ai,  1991;  Rosakis,  1993)  described 
below. 


7.1  Experimental  technique  (transmission  CGS) 

Consider  a  planar  wavefront  normally  incident  on  an  optically  and  mechanically  isotropic, 
transparent  plate  of  initial  uniform  thickness  h  and  refractive  index  n.  As  shown  in 
Figure  7,  the  specimen  occupies  the  (xj ,  x2)  plane  in  the  undeformed  configuration.  When 
the  specimen  undergoes  any  kind  of  deformation  (static  or  dynamic),  the  transmitted 
wavefront  can  be  expressed  as  S(xi,x2,x3)  =  X3  +  AS(xj,x2)  =  constant,  where  AS  is 
the  optical  path  change  acquired  during  refraction.  As  discussed  in  detail  by  ROSAKIS 
(1993),  AS  is  related  to  the  deformation  state  by  the  relation, 

AS(xj,x2)  =  2h(n  -  1)  /  C33d(x3//i)  -f  2h  f  And(x3/h)  .  (7.1) 

Jo  Jo 


The  first  term  of  equation  (7.1)  represents  the  net  optical  path  difference  due  to  the 
plate  thickness  change  caused  by  the  strain  component  C33.  The  second  term  is  due  to 
the  stress  induced  change  of  refractive  index  of  the  material.  This  change  in  refrative 
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index  An  is  given  by  the  Maxwell  relation, 


An  =  D\  (an  -f  <^22  4-  <733)  , 


(7.2) 


where  D\  is  the  stress  optic  coefficient  and  a \j  are  Cartesian  components  of  the  nominal 
stress  tensor.  The  above  relation  is  strictly  true  for  isotropic,  linearly  elastic  solids.  For 
such  solids,  the  strain  component  633  can  also  be  related  to  the  stresses,  and  equation 
(7.1)  then  becomes: 


AS(xx,x2)  =  2hca  '  j(c„+<722)  \-D2  }  (*3/*)  ,  (7.3) 


where 


ca  =  D{ 


Hn~  1) 


(<7n  4-  <r2  2)y 

u{n  —  1) 

vDi  4-  — -7; - 

D2  - - E- 
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Dx 


and  £,  1/,  and  c„  are  the  Young’s  modulus,  Poisson’s  ratio  and  stress  optical  coefficient 
of  the  material,  respectively. 


A  schematic  of  the  experimental  apparatus  is  also  shown  in  Figure  7.  VVhen  the 
transmitted  wavefront  emerges  from  the  specimen  after  being  distorted,  it  passes  through 
two  high  density  gratings,  G\  and  G2  of  pitch  p,  separated  by  a  distance  A.  The  gratings 
have  their  rulings  parallel  to  either  the  Xj  or  x2  directions.  The  action  of  the  gratings 
is  to  displace  (shear)  the  diffracted  beam  and  recombined  it  with  itself,  thus  creating 
an  interferogram  after  G2.  The  filtering  lens  L  processes  the  light  emerging  from  G2 
and  its  frequency  content  (diffraction  spots)  is  displayed  on  the  back  focal  plane  of  L. 
By  physically  blocking  all  diffraction  orders  except  for  either  the  ±1  orders,  information 
regarding  the  gradient  components  of  A5(xj,  x2)  along  either  the  xi  or  x2-axis  is  obtained 
on  the  image  plane.  The  camera  is  kept  focused  on  the  specimen  plane.  For  grating  rulings 
perpendicular  to  the  xa-axis,  the  resulting  fringe  pattern  is  proportional  to  d(A S)/dxa, 
{1,2}. 

A  first  order  analysis  described  by  TlPPUR  et  al.  (1991),  or  a  higher  order  Fourier 
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Invariably  a  near-tip  three-dimensional  region  will  exist  in  any  real  specimen  geometry. 
However,  outside  this  three-dimensional  zone,  a  plane  stress  approximation  will  be  valid. 
A  numerical  study  of  each  particular  specimen  configuration  is  needed  to  identify  the 
extent  and  exact  location  of  such  a  plane  stress  region.  Such  a  calculation  has  been 
performed  by  LEE  and  ROSAKIS  (1992)  for  a  three  point  bend  bimaterial  specimen.  A 
rather  large  two-dimensional  plane  stress  region  was  seen  over  a  significant  portion  of  the 
specimen.  In  this  region,  a33/r'(orn  -f  <r2 2)  (a  measure  of  three-dimensionality;  tends  to 
zero.  For  points  ouside  the  three-dimensional  region  n  +  022)  0)?  the  optical 

path  difference  in  equation  (7.3)  will  simplify  to 

AS(xux2)  ~  cvh{au(xux2)  +  022(21, *2)}  ,  (7.5) 

where  da  and  d2 2  are  thickness  averages  of  the  stress  components  in  the  plate. 


As  a  result,  for  points  outside  the  near-tip  three-dimensional  region,  the  CGS  patterns 
assume  a  simple  interpretation  in  terms  of  two-dimensional  stress  field  app>  ximations. 
In  particular,  equations  (7.4)  and  (7.5)  now  indicate  that  fringes  obtained  from  regions 
surrounding  the  three-dimensional  zone  can  be  related  to  the  in-plane  gradients  of  da  +b22 


as  follows: 


,d(<Tn+v 22)  mP  ,d(<rn+cr22)  np  _  ,  cx 

Ceh-^-r - 1  =  —  ,  cah - - - =  -T-  ,  m,  n  =  0,  ±1,  ±2,  •  •  •  (7.b) 

dij  A  dx2  A 

where  in  the  case  of  transmission  c„  is  the  stress  optical  coefficient  of  the  material  (e.g. 


PMMA). 


51 


7.2  Experimental  set-up  and  procedure 


Bimaterial  specimens  used  in  the  dynamic  experiments  are  of  the  three  point  or  one 
point  bend  configuration  and  are  made  from  9mm  thickness  sheets  of  commercially  avail¬ 
able  poly-methylmethacrylate  (PMMA)  (material-1)  and  A1SI  4340  steel  (material-2). 
The  bonding  procedure  is  outlined  in  TlPPUR  and  RosaKIS  (1990).  A  bond  strength 
calibration  experiment  was  also  performed  in  that  study,  demonstrating  that  the  bond 
toughness  was  at  least  as  much  as  that  of  a  homogeneous  PMMA  specimen.  This  fact 
testifies  to  the  strength  of  the  bond  and  becomes  important  in  the  discussion  of  the 
dynamic  experiments  presented  bellow. 

The  bimaterial  specimens  have  either  a  pre-cut  edge  notch,  or  a  sharp  pre-crack  of 
length  25mm  along  the  interface.  The  specimens  are  either  impact  loaded  in  a  drop 
weight  tower  (Dynatup-8100A)  or  a  high  speed  gas  gun.  After  the  impact  event,  the 
crack  propagates  dynamically  along  the  interface.  The  transmission  CGS  technique  in 
conjunction  with  high  speed  photography  is  used  to  record  dynamic  fields  around  the 
crack-tip  (only  on  the  PMMA  side,  of  course).  A  rotating  mirror  high  speed  camera 
(Cordin  model  330A)  is  used.  A  Spectra- Physics  Argon-ion  pulse  laser  (model  166)  is 
used  as  the  light  source.  By  using  short  pulses  of  30nsec  duration,  we  are  able  to  freeze 
even  the  fastest  of  running  cracks  and  thus  produce  a  sharp  interference  pattern  during 
crack  growth.  The  interframe  time  (controlled  by  the  interval  between  pulses)  is  typically 
set  at  l^sec  for  a  total  recording  time  of  80/xsec.  The  laser  pulsing  is  triggered  by  a  strain 
gauge  on  the  specimen  that  senses  the  impact. 

True  symmetric  one  or  three  point  bend  loading  cannot  be  achieved  since  it  is  ex¬ 
tremely  difficult  to  apply  the  impact  load  exactly  on  the  interface,  which  is  very  thin. 
In  addition  since  the  wave  speeds  of  PMMA  and  steel  are  vastly  different,  the  loading 
history  at  the  crack-tip  would  be  completely  different  if  the  specimen  were  impacted  on 


52 


the  PMMA  or  the  steel  side.  Thus  it  was  chosen  to  impact  the  specimen  a  small  distance 
(7mm)  into  the  steel  side  of  the  bond. 

A  sequence  of  high  speed  interferograms  from  a  PMMA/steel  test  is  shown  in  Figure  8. 
This  is  a  three  point  bend  test  conducted  in  a  drop  weight  tower.  The  impact  speed  was 
4m/sec.  When  the  crack  initiates  (t  =  0/isec),  intense  stress  waves  emanate  from  the 
crack-tip.  These  w’aves  are  visible  in  Figure  8  as  discrete  kinks  in  otherwise  smooth  fringes 
and  as  circular  lines  centered  at  points  along  the  crack  line  (see  frames  at  t  =  16.5/isec 
and  t  =  23/isec).  This  observation  is  a  reliable  sign  of  a  highly  dynamic  event,  as  will  be 
discussed  later. 


7.3  Analysis  of  experimental  data 

In  subsequent  sections  we  shall  present  an  analysis  of  CGS  interferograms  of  dynamic 
bimaterial  specimens  first  using  a  /^-dominant  assumption  and  then  using  the  higher 
order  transient  field  described  in  section  4- 

7.3.1  Singular  field  (A'd-dominance) 

The  governing  relations  for  CGS  (7.6)  can  be  used  to  estimate  fracture  parameters  from 
points  outside  the  three-dimensional  zone  of  a  given  interferogram.  One  could  expect 
that  the  plane  stress  region  surrounding  the  near-tip  three-dimensional  region  would  be 
well  described  by  the  most  singular  term  in  the  asymptotic  expansion  for  stress,  i.e.  that 
a  /^-dominant  region  would  exist  somewhere  around  the  crack-tip.  This  is  something 
to  be  verified  though  and  should  not  be  taken  for  granted,  especially  in  regions  relatively 
far  from  the  crack-tip  or  in  experiments  showing  transient  effects  (e.g.  rapidly  changing 
crack-tip  speed).  In  such  cases  the  deformation  field  around  the  crack-tip  may  be  better 
described  by  a  higher  order  analysis. 
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As  was  stated  earlier  (see  section  4),  for  cracks  propagating  dynamically  under  steady 
state  conditions  in  bimaterial  specimens,  Yang  et  al.  (1991)  and  the  first  part  of  the 
present  analysis  observed  that  near  the  crack-tip  the  stress  field  assumes  the  form, 


where  (r,  0)  are  polar  coordinates  of  a  coordinate  system  translating  with  the  crack-tip  at 
speed  v,  and  Kd  is  the  complex  dynamic  stress  intensity  factor.  The  material  mismatch 
parameter  e  =  e(v)  is  now  a  function  of  crack-tip  speed  and  of  the  elastic  moduli  of  the 
materials  of  the  bimaterial  system.  Analytical  expressions  for  and  dj^  are  given  by 
Yang  et  al.  (1991). 

By  using  equation  (7.7)  and  after  some  algebraic  manipulations,  dn  -f  a  22  can  be 
written  as 


au+a22  =  -^L|(l-fa^-27/a,)et(,r~Si)cos^-$(t)-elnr/^ 

+  (l  +  q]  +  217a,)  cos  +  $(*)  +  £  In  r^j  j 


(7.8) 


where 


(a?  -  a?) 

A"(0 

D(v)  cos 

b(£7T) 

,  Kd(t)  =  Ki(t)  +  ,  *(()  =  tan'1  , 


and  at,,,  r/t„  and  0{}3  have  been  defined  in  previous  sections.  The  mismatch  parameters 
r)  and  c  are  functions  of  crack- tip  speed  and  of  material  properties.  These  functions  are 
given  in  section  5  and  appear  in  Figures  2  and  5,  respectively.  Note  that  equation  (7.8) 
is  the  first  part  of  equation  (4.6)  in  section  4.  The  field  quantity  of  interest  in  analyzing 
the  CGS  patterns  for  material-1  is  cahd(ou  -f-^22)/ dx\.  By  differentiating  equation  (7.8) 
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with  respect  to  X\,  we  have 


d(an  +  £22)  =  Crhr^e'^-^Ajt) 

dxi  2>/2tt 

x  |  -  (l  +  a]  -  2 r?as)  e2e(*~e,)  cos  -  $(f )  -  e  In 


-  (l  +  a]  +  2»7a3)  cos  +  $(*)  +  clnr/j 

+2e  (l  +  a2  —  2r;a,)  e2d*-*i)  sin  —  $(t)  -  t  In  r^j 

-2c  (l  +  a]  +  2t?qs)  sin  +  $(t)  +  e  In  r,  j  j 
where  A(t)  is  as  defined  in  equation  (7.8)  and  0  <  Oi  <  it. 


)  , 


(7.9) 


From  the  above  discussion  it  becomes  obvious  that  extraction  of  parameters  like  Kd 
is  now  possible  provided  that  experimental  data  are  gathered  from  a  region  near  the 
moving  crack-tip  characterized  by  the  structure  presented  in  equations  (7.8)  and  (7.9). 
In  a  laboratory  specimen  of  finite  size  where  transient  effects  may  be  important,  the 
field  may  not  be  /^-dominant  and  the  use  of  a  higher  order  analysis  may  be  necessary. 
The  necessity  of  a  higher  order  analysis  in  the  interpetation  of  optical  data  from  crack 
growth  in  homogeneous  specimens  was  demonstrated  by  FREUND  and  Rosakis  (1992) 
and  KRISHNASWAMY  and  ROSAKIS  (1991).  An  equivalent  analysis  for  a  transiently 
propagating  interfacial  crack  has  been  provided  in  previous  sections  and  its  effect  on 
data  interpretation  is  discussed  in  the  next  section. 


7.3.2  Higher  order  transient  analysis 

In  section  4,  a  higher  order  expansion  for  the  trace  of  the  stress  tensor  in  plane  stress 
is  shown  in  equation  (4.6).  By  differentiating  with  respect  to  the  xx  coordinate,  we 
obtain  a  relation  for  the  xj-gradient  of  <7n  +  <722,  which  is  relevant  to  the  analysis  of  CGS 
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interferograms, 

iV/22^1  =  | AW!  {noWcos(C In n)+no(^)sin(e In r,)}rr 

2  fi(af-a])  l  i 


4-  e  |n,i(5/)cos(clnr;)4-  (0/)  sin(elnr,)|  r,-1/2(ln  r/)J 

4-  cos(c  In  r/)+  fh  (0/)  sin(e  In  r;)|  r^l/2  In  r, 

4-  |n„(0|)  cos(e  In  r,)+  jjtt  (0/)  sin(e  In  r/)|  r,~1/2 
+  |A(0I  |n2(^)cos(clnr/)+  fl2  (0/)  sin(c In  r,)j  rf1/2 


,  (7.10) 


+  0(rt) 


where 


no(0,)  =  a0(t)e-tff>  cos  ^  -  *(*,))  +  M*V*  cos  ^  +  *(  A))  . 
no  (0;)  =  ao(t)e~te‘  sin  -  $(A) j  -  bo(t)ete‘  sin  4-  $(  A)  j  , 
n2(0,)  =  a0(Oe-te,cos^|-$(A))+MOceff'cos^|  +  $(A)j, 
n2  (0/)  =  a0(t)e-t6'  sin  -  $(A)  j  -  M*K®‘  sin  4-  $(A)^  , 

H»(0,)  m  |i2d(0k~te,cos^|-<>(e^-|ea0|c<Sicos^  +  $(G^, 

L  (0/)  =  |Qd(«)|e"<9'  sin  -  $(<2d)j  +  |  2d  {t)\e(9t  sin  4-  $(Cd))  , 
nt(0,)  =  |A(0|c“rf' cos  (f  ~  *(A )j  ~  I  At  {t)\et6'  cos  ^  -  $(A)j 

4-  |Bt(0|e-£<’'  cos  -  I  Bt  (01^®'  cos  “  *  A)) 

4-  2c  j  jQd(0ie*"e*‘ sin  “  *(&))  -  !  Qd  (0kt<>(  sin  "  *{QdYj  |  9, 
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fit  (ft)  =  IA(0|e-£fl'  Sin  -  *(.**))  +  I  At  (0|e^‘  sin  -  d>(A) j 
+  \Bt(t) \e~ie‘  sin  -  $(Bt)J  +  |  Bt  (0 |c'a'  sin  -  $(Bt)  j 

-  2e  jlG^le-*  cos  -  4>«2„))  +  \  Qd  (t)\e<*‘  cos  -  *(Qd)J  }  ft, 

n«(ft)  =  |At(^)|e_t5'  COS  -  |  Att  (0|etfii  cos  ^y  -  ${Att)^ 

+  |B«(*)|e-'®'  cos  ^ - y  -  $(#«)^  -  |  Bu  (t) \et6‘  cos  ^^y  -  $(Btt)j 
+  I Cu(t)}e~c61  cos  ^^y  -  $(Ctt)j  -  |  Ctt  (<)|eeff|  cos  ^y  -  $(C«)^ 

+  |l-4t(0le  <fi<  s*n  ~  $(A)j  -  |  At  (f)|ec®'sin  -  $(.4,)^  j  ft 

+  ||S1(<)|c-t9‘  sin  -  $(Bt)  j  -  |  Bt  {t)\e‘B‘  sin  -  *(Bt)  j  |  9, 

-  i  ||(2d(0|c-tfl'  cos  -  *(&))  -  |  Qi  (<)!*“'  cos  J  6f, 

flu  (ft)  =  \Mi)\e-t6'  sin  -  *M„))  +  I  Att  (t)\et6‘  sin  -  $U„)) 

+  |B«(t)|e-^  sin  ^  -  *(B«i))  +  |  Bu  (*) Je*  sin  ^  -  *($«)) 

+  |C«(0|e'£<?'  sin  (—■  -  $(C«)j  +  |  Ctt  (0 |e£®'  sin  ^y  -  $(£«)  j 

-  |lA(<)|e“‘9'  cos  ^y  -  +  |  At  (<)ketft  cos  ^y  +  $M<)^  |  ft 

~  j|Bi(<)|e“rt|  cos  ^y  -  $(Bt)j  +  |  Bt  (t) \e<9‘  cos  -  $(Bt)j  J  ft 

-  c  jl&MI^'sin  ^y  -$(£>«,) j  +  |  Qd  (t)\et0t  sin  ^y  -  $(Qrf)j|ft2, 
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The  functions  of  time  Ao{t),  A2{t),  Au(t),  •  •  •,  that  appear  in  the  above  expressions,  are 
related  to  functions  Ao(t)i  -''MO*  Att(t),  ■  •  in  equation  (4.5)  by 

Ao (0  =  (~ 2  *  *e)  (*)  ’  ^(O  =  (2  +  *e)  ^  ' 

Au{t)  =  +  At(t)  +  Btt(t)  , 

An  (t)  =  (|  -  ie)  A«  (0+  A,  (0+  £»  (<)  , 


B«[t)  = 

(-!  +  “) 

1  Btt(t)  +  Bt(i)  +  2 Cu(t)  , 

= 

(~H 

1  Bu  (*)+  Bt  ( t )  +  2  Ctt  ( t )  , 

C«(t)  = 

,  •  /  3  \  • 

(-2+’£) 

1  c»(0  1  Ctt  (t)  —  y--  —  itj  Ctt  (t)  1 

—  (2  *e)  ^«(0  +  2cfi,*(t)  4-  Bt(t)  , 

A  (0  =  (|  -  «)  Att  (t)  +  2c  nd  (0+  Bt  (t)  , 
fit(0  =  +  «)  *«(<)  >  Bt  (0  =  (~  -  «)  Bt  (t)  , 

Qd(0  =  (2  l€)  (0  =  (<J  “  *c)  fk  (0  • 

This  gradient  contains  4  orders  in  rj.  They  are  rj^2 ,  r/-1^2(lnr/)2,  rj_I^2lnr/,  and  r,-1^2. 
It  also  contains  28  undetermined  constants.  The  first  two  constants  |,4o|  and  4>(>to)  are 
related  to  \Kd\  and  $  (or  Kd,  Kd)  of  the  expression  of  Yang  et  al.  (1991)  (see  equation 
(7.9)).  In  fact  the  most  singular  term  of  equation  (7.10)  reduces  to  equation  (7.9).  Under 
steady  state  conditions,  equation  (7.10)  reduces  to  an  expression  with  4  terms  which 
are  identical  to  the  first  4  terms  of  the  higher  order  steady  state  expression  derived  by 
DENG  (1992).  The  transient  contributions  to  the  expression  for  the  gradient  (7.9)  are 
those  that  exbihit  an  r,-1^2(lnrj)2  and  rj"1^2  In  rj  radial  dependence.  It  is  worth  noting 
that  most  of  these  transient  terms  are  multiplied  by  the  quantity  e,  the  rate  of  change  of 
the  ocsillatory  index  with  time  (c  =  e'(v)v)-  Thus,  to  a  certain  extent,  e  is  a  measure  of 


transience  of  the  propagating  crack.  If  k  —  0,  most,  but  not  all  transient  terms  disappear. 
Those  that  remain  are  those  related  to  the  rate  of  change  of  the  complex  stress  intensity 
factor.  Note  that  it  is  possible  for  e  to  be  small  even  if  a  large  acceleration  exists,  but 
e'(v)  is  small.  Conversely  it  is  possible  to  have  a  large  e  corresponding  to  small  v  but 
large  e'(u).  It  should  be  noted  that  d(v)  tends  to  infinity  as  v  tends  to  C/j*,  see  Figure  5. 
Whether  or  not  c  can  be  used  as  a  reliable  measure  of  transience  will  be  investigated  in 
the  subsequent  section. 

It  is  clear  at  this  point  that  analysis  of  the  fringe  patterns  obtained  from  a  dynamic 
experiment  can  be  made  using  either  equation  (7.9)  or  equaton  (7.10).  The  choice  of  one 
or  the  other  depends  on  whether  a  region  of  /^-dominance  has  been  estalished  somewhere 
outside  the  near-tip  three-dimensional  zone.  Use  of  either  equation  allows  estimation  of 
the  time  variation  of  the  relevant  parameters.  This  is  done  by  performing  a  least  squares 
fitting  procedure  to  data  points  digitized  from  the  CGS  interferograms  obtained  during 
an  experiment.  Of  course  the  crack-tip  speed  n(t)  is  measured  independently.  There  are 
2  undetermined  parameters  in  equation  (7.9)  and  28  undetermined  constants  in  equation 
(7.10). 

7.4  Results  and  discussion 

The  velocity  and  acceleration  histories  corresponding  to  the  sequence  of  photographs 
in  Figure  8  are  shown  in  Figures  9(a)  and  9(b).  This  is  a  test  performed  in  a  drop 
weight  tower  under  the  relatively  small  impact  speed  of  4m/sec.  Indeed  the  terminal 
speed  in  this  test  seem  to  be  about  90%  of  the  Rayleigh  wave  speed  of  PMMA,  see 
Figure  9(a).  In  contrast,  previous  experience  with  dynamic  crack  growth  in  homogeneous 
PMMA  specimens  of  the  same  configuration  show  a  maximum  speed  of  about  0.35c/j*. 
Note  also  that  in  this  particular  bimaterial  case  there  is  a  very  large  crack-tip  acceleration 
(approximately  107g,  where  g  is  the  acceleration  of  gravity)  immediately  after  the  crack 
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initiates,  see  Figure  9(b).  This  would  suggest  that  transient  effects  would  be  present  close 
to  initiation  ( t  =  0//sec).  As  was  mentioned  earlier  the  rate  of  change  of  the  oscillatory 
index  with  time  (e)  may  be  considered  a  partial  measure  of  transience.  For  the  same 
test  as  Figure  8,  we  have  plotted  t  and  i  versus  time  in  Figures  10(a)  and  10(b).  In 
Figure  10(b),  i  exhibits  a  local  maximum  at  about  t  =  10/isec  after  initiation.  It  then 
starts  increasing  again  after  25/zsec.  At  short  times  after  initiation,  e'(v)  is  close  to 
zero  although  v  is  large  (107g).  This  accounts  for  the  initially  low  values  of  e.  In  this 
regime  transient  effects  are  demonstrated  through  large  changes  in  the  complex  dynamic 
stress  intensity  factor.  As  time  increases  the  combinaton  of  e'(v)  and  v  results  in  a  local 
maximum  in  L  At  later  times  ( t  >  25//sec)  and  as  the  crack-tip  velocity  approaches  the 
Rayleigh  wave  speed  of  PMMA,  i  increases  again. 

To  demonstrate  the  need  of  a  transient  analysis  in  interpreting  experimental  data, 
let  us  now  attempt  to  analyze  the  frame  of  Figure  8  at  t  =  9.5psec.  This  corresponds 
to  a  local  maximum  value  of  i  in  this  particular  test.  By  following  the  fitting  proce¬ 
dure  described  in  section  7.3.2,  we  can  obtain  the  coefficients  of  either  equation  (7.9)  or 
equation  (7.10).  The  result  of  such  a  fit  for  the  /^-dominant  field  (equation  (7.9))  is 
shown  in  Figure  11(a).  The  diamonds  are  digitized  data  points  from  the  interferogram 
at  t  =  9.5psec.  The  solid  line  is  the  contour  of  the  quantity  d(bn  +  022)!  dx\  calculated 
numerically  by  using  the  results  for  Kd  from  the  fit  generated  by  the  same  data  points. 
As  can  be  clearly  seen,  equation  (7.9)  cannot  represent  the  data  to  any  reasonable  extent. 
The  deformation  field  of  this  particular  picture  therefore  is  nowhere  near  /^-dominant. 
In  fact  the  main  feature  which  is  that  the  fringes  vertically  approach  the  interface  cannot 
be  captured  at  all  by  equation  (7.9).  The  result  of  the  fit  of  the  transient  higher  order 
field  (equation  (7.10))  derived  earlier  is  shown  in  Figure  11(b).  The  data  points  are  ex¬ 
actly  the  same  as  before  and  the  solid  line  is  the  result  of  the  fit.  Clearly  the  fit  is  very 
good  over  a  large  area  of  the  specimen.  All  features  of  the  field  are  successfully  captured 
by  equation  (7.10).  This  shows  that  the  A'rf-dominant  analysis  cannot  be  used  for  cases 
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where  e  is  high. 


To  further  investigate  the  effect  of  i  on  the  interpretation  of  optical  data,  we  chose  to 
analyze  an  interferogram  corresponding  to  the  minimum  value  of  e  within  the  duration 
of  the  test.  This  occurs  at  t  =  23/zsec.  Figure  12(a)  shows  the  result  of  the  /^-dominant 
fit  to  the  experimental  data.  As  the  crack-tip  is  approached,  equation  (7.9)  seems  to 
adequately  describe  the  experimental  measurement.  However,  as  the  distance  from  the 
crack-tip  is  increased,  Kd- dominance  is  lost.  Nevertheless,  the  lack  of  A'rf-dominance  in 
Figure  12(a)  (e  ~  1.0 x  102sec-1)  is  not  as  dramatic  as  in  Figure  11(a)  (e  ~  1.2  x  lO^sec-1). 
Figure  12(b)  shows  the  result  of  the  fit  of  the  transient  higher  order  field  to  the  same 
experimental  data  as  Figure  12(a).  The  fit  is  now  much  better  over  the  whole  range 
of  radii.  The  above  observations  show  that  in  general  a  transient  analysis  of  data  is 
necessary  if  fracture  paxameters  such  as  Kd  are  to  be  obtained  with  confidence. 

7.5  Transonic  terminal  speeds 

The  next  cycle  of  experimentation  involved  bimaterial  specimens  loaded  at  higher  loading 
rates  than  in  a  drop  weight  tower.  This  was  achieved  by  using  a  high  speed  gas  gun.  A 
one  point  bend  impact  geometry  was  used.  Again  the  issues  of  crack-tip  loading  history, 
as  dependent  upon  PMMA  or  steel  side  impact,  arise.  It  was  chosen  to  impact  the 
specimens  on  the  steel  side,  to  remain  consistent  with  the  drop  weight  tower  tests.  The 
gas  gun  projectile  was  50mm  in  diameter  and  the  impact  velocity  was  20m/sec,  thus 
resulting  in  considerably  larger  near-tip  loading  rates  than  in  the  drop  weight  device.  A 
sequence  of  interferograms  from  such  a  test  is  shown  in  Figure  13.  Its  corresponding  t>(<), 
v(t),  c(f)  and  e(t)  plots  are  shown  in  Figures  14(a),  (b),  and  15(a),  (b).  In  general  terms 
the  results  are  similar  to  those  obtained  from  the  drop  weight  tower  experiments.  A  main 
difference  is  that  the  speed  and  acceleration  are  much  higher.  In  fact  the  crack-tip  speed 
seems  to  exceed  the  Rayleigh  wave  speed  of  PMMA  after  a  relatively  short  time.  In  some 


61 


cases  (as  in  Figure  14)  the  velocity  even  exceeds  the  shear  wave  speed  and  approaches 
the  longitudinal  wave  speed  of  PMMA,  thus  entering  the  transonic  speed  range  for  the 
PMMA  side. 

For  a  crack  speed  less  than  the  Rayleigh  wave  speed,  we  can  repeat  a  fitting  procedure 
exactly  as  before.  For  the  frame  at  t  =  8/xsec  in  Figure  13,  the  result  of  such  fit  is  shown 
in  Figure  16.  Here  the  white  lines,  obtained  from  plotting  the  field  of  equation  (7.10) 
using  the  values  of  the  fitted  parameters,  are  superposed  on  the  actual  picture  (instead 
of  the  digitized  points  as  in  Figures  11  and  12).  The  illustration  is  the  same  though,  i.e. 
that  a  transient  field  is  necessary  to  describe  a  picture  such  as  this  which  corresponds  to 
a  high  e  and  acceleration. 

Unfortunately  given  the  existing  theoretical  analyses,  we  do  not  have  the  tools  to  fit 
any  field  to  interferograms  having  a  speed  in  the  transonic  range  for  PMMA  (c^  <  v  < 
c|^).  These  large  speeus  were  observed  in  a  number  of  tests  involving  one  point  bend 
interfacial  specimens  containing  sharp  pre-cracks  lying  along  the  interface.  When  a  spec¬ 
imen  containing  a  blunt  starter  notch  was  impacted,  recorded  crack-tip  terminal  speeds 
were  even  higher;  at  some  cases  approaching  the  longitudinal  wave  speed  of  PMMA.  Such 
a  velocity  history  is  given  in  Figure  17.  Here  the  maximum  crack-tip  speed  is  estimated  to 
be  O^cJ1^.  These  observations  are  very  interesting  because  to  our  knowledge  no  evidence 
of  transonic  or  supersonic  crack  propagation  has  ever  been  seen  in  homogeneous  materi¬ 
als  even  though  a  large  number  of  the  theoretical  studies  exist  on  the  subject  (FREUND, 
1990).  It  is  believed  that  transonic  crack  growth  is  possible  in  a  bimaterial  situation 
because  of  an  energy  transfer  mechanism  from  the  stiffer  to  the  softer  material.  It  can  be 
seen  in  Figure  13  that  the  nature  of  the  fringes  changes,  approximately  around  the  time 
at  which  the  crack-tip  speed  exceeds  the  Rayleigh  wave  speed.  A  sequence  corresponding 
to  the  same  test  whose  velocity  is  shown  in  Figure  17  (blunt  starter  notch)  is  presented 
in  Figure  18.  In  these  pictures,  we  see  an  even  more  drastic  change  in  the  nature  of  the 
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fringe  patterns  as  the  crack-tip  speed  exceeds  both  Rayleigh  and  shear  wave  speeds.  To 
see  this  effect  clearly,  compare  the  second  frame  in  Figure  18  to  the  sixth  frame.  Finally, 
additional  visual  proof  of  the  existence  of  large  transient  effects  is  shown  in  Figure  19. 
We  are  now  in  the  process  of  developing  an  analysis  for  the  propagation  of  an  interfacial 
crack  at  speeds  exceeding  c^.  It  is  hoped  to  be  able  to  predict  fringe  patterns  as  those 
observed  in  Figures  13,  18,  and  19. 

8  Conclusions 


Experimental  observations  of  high  speed  (transonic  terminal  speeds)  and  high  accelera¬ 
tion  (108m/sec2)  crack  growth  in  PMM A/steel  interfaces  are  reported  for  the  first  time. 
Motivated  by  these  observations,  a  fully  transient  higher  order  asymptotic  analysis  of 
dynamic  interfacial  crack  growth  is  performed.  This  analysis  is  valid  for  crack-tip  speed 
in  the  range  0  <  v  <  (c^  is  the  shear  wave  speed  of  PMMA).  Explicit  expressions 

for  stresses  are  provided.  In  addition  to  the  classical  r-1/2,  r°,  and  r1/2,  •  -  •,  terms  of 
steady  state  expansion  for  the  stresses,  new  transient  contributions  of  order  r1^2  In  r  and 
r1/2(ln  r)2  appear.  The  structure  of  the  near-tip  field  obtained  by  the  analysis  is  found 
to  describe  well  the  experimentally  obtained  stress  fields.  For  subsonic  crack  growth,  the 
experiments  demonstrate  the  necessity  of  employing  the  fully  transient  expression  in  the 
analysis  of  optical  experimental  data.  Terminal  speeds  of  up  to  90%  of  the  plane  stress 
dilatational  wave  speeds  of  PMMA  are  observed. 
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Appendix 


by 


I:  Definitions  and  properties  of  matrices  used  in  section  3 

Let  JPfe,  Qfc,  U fc,  and  be  defined  as  in  page  12,  section  3,  and  L k  and  Lk  be  given 

vkQ:1 


H=Li  -L2  . 


Lk  =  UkPkx  , 

* 

Matrices  H  and  H  are  defined  as, 

H  =  Li—  L2  , 

By  algebraic  calculations,  it  can  be  shown  that  for  k  6  {1,2}, 
Lk  = 

where 


(^ll)fc  (^12)  k 

* 

'  (/«)*  -(/«)*  ' 

.  (^21  )fc  (hl)k  . 

,  Lk  = 

.  “  (All)*  (All)*  . 

f2a,a4  -  (1  +  a2s)\  ,,  *  fa,(l-a2)]  [^(l-a^)} 

— im — h’  (W*  -  r^n  ■ 


and 


Therefore, 


D{v)  =  4aja,  -  (1  +  a2)2  . 


JFf  = 


h\2 
h2\  h\\ 


hn  —hi2 
—  h2\  h\\ 


where 

hn  =  (hi)i  ~  (hih  •>  hn  =  (^12)1  *f  (*12)2  »  ^21  =  (*21)1  + 

Notice  that 

H  H-H  H  =  (ft2!  -  huh2l)I , 

where  /  is  the  2x2  identity  matrix.  Thus, 


1  *  *  1 
if1  =  H  ,  H  =  : — H 


^11  ~  h\2h2i 


^11  —  ^12^21 


69 


Also,  it  can  be  shown  that 


Lk  Lk=Lk  Lk  ,  k  (E  {1,2}. 

A  sequence  of  operator  definitions  follows.  These  are  related  to  the  analysis  in  section 
3.2.  Let  p(t)  and  q(t)  be  two  real  functions  of  time  t  and  define  the  vector  operators 

dk{p(t),q(i)}  =  {  Ddpit)},  D,{q(t)}  }Tk  , 

kk{p(t),q(t)}  =  {  K,(t)p(t),  Ka(t)q(t)  , 

^(«)>  =  {  Bi(t)p(t ),  Ba{t)q(t)  )l  , 

where  operators  £>/,,{■}  and  functions  Kt,a{t)  and  Bi,a(t )  have  been  defined  in  section  3.2. 
With  the  above  definitions, 

tk{p{t),  q{t)}  =  (3  +  2 it)  ( LkMk  -  I)  dk{p(t),  q(t)}  +  2e  { LkMk  -  I)  **{p(<),  ?(*)} 

+  2  [(1  -f  2 it)  ( LkMk  -I)  +  LkPk-  £*]  bk{p(t),  q(t)}  , 

*  * 

where  Af*.,  Pk,  and  Uk  have  also  been  defined  in  section  3.2.  In  addition,  for  any  given 
operator 

*Mp(*),9(*)}  =  {  ?(*)},  mfc2){p(0,?(0}  }  , 

the  associated  operator  m*  {?(*),?(*)}  is  defined  as 

{p(t),<i{t)}  =  { "4l){p(<)> $(*)}>  -TOLa){fK0.9(0}  }T  • 

Also  vectors  7,  and  c,  are  defined  as 
0  =  ti{a0(f),co(t)}-  t2  {io(t),do(<)} 

7  =  *1  {M*Mo(0} 

£  =  (3  +  2te)  [(LjMi  -  /) *, {ao(t), <*(*)}  -  (12  M 2  -  /)  fc2  {6o(t), 4,(0}] 

<  =  (3  +  2ie)  [(L2M2  -  /)*2{«o(*W<)}  -  Afi  -  /)  *1  {M<Mo(0}] 
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and  operators  Wdk{p(t),q(t)}  an<^  tDtjt{p(<),^(<)}  as 

Wdk{p(t), q(t)}  =  (|  +  *e)  [2 PlxMk  +  (|  +  te)  /]  **r{p(<),  <?(*)} 

*><*  {?(*)>?(<)}  =  (|  +  «)  +  iej  /]  dk{p(t),q{t)} 

+  2e  [P^Mt  +  (1  +  ie)l]  kk{p(t),q(t)} 

+  [2 P;1  Pk  +2(1  +  2 it)P~kxMk  -  Q  +  e2)  /]  bk{p(t),q(t)} 

II:  Solution  of  the  Riemann-Hilbert  problem 

Consider  the  problem  formulated  as  following:  Find  a  function 

«W  =  ( #,M,  «,(*)  )T  , 

z  —  7/1  +  *772,  which  is  analytic  in  the  whole  z-plane  except  along  the  branch  cut  —  00  < 
Vi  <  0,  r}2  =  0)  and  satisfies  the  equation, 

H «+(,,)-  JT«-(ih)  =  <c(„)  ,  V„<0,  (1) 

where  JEf  and  JFf  are  2x2  matrices,  defined  in  Appendix  I,  and 

*(>71)  =  (  *2(771)  )T  , 

with  K\  and  k2  are  known  functions  of  r)\.  Near  the  origin,  function  0(z)  should  satisfy 
the  requirement  that 

1  e(z)  I  =  o ( |*r ) .  »W-0,  (2) 

for  some  real  number  q,  and  generally,  a  >  —  1. 

In  order  to  obtain  the  solution  to  the  above  Riemann-Hilbert  problem,  the  eigen- values 

m 

and  eigen-vectors  of  H ,  and  H  need  to  be  studied  first.  By  solving  the  equation 

det  {H  -  A/}  =  0  ,  (3) 
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where  I  is  the  identity  matrix,  the  eigen-values  for  H  are  found  to  be 


^1,2  —  ^11  ±  \Jh12h21  ■  (4) 

The  expressions  of  functions  hn,hx2,  and  h2 1  in  our  problem  are  dependent  upon  the 
mechanical  properties  of  the  constituents  of  the  bimaterial  system  and  the  speed  of  prop¬ 
agation  of  the  interfacial  crack.  hn,h12,  and  h2i  ensure  that  the  eigen-values,  A]  and  A2, 
are  real,  provided  that  the  crack-tip  speed  is  less  than  the  lower  Rayleigh  wave  speed  of 
the  bimaterial.  The  corresponding  eigen-vectors  are 

™(1’2)  =  (  1,  ±T})J  ,  (5) 

where  the  parameter  77  is  defined  by 


It  can  be  shown  that  the  eigen- values  for  H  are  the  same  as  those  for  H,  which  are  given 
in  (4),  while  the  corresponding  eigen- vectors  are 


*  /  1  \T 

™  =  (  1,  ) 


Define  the  matrix  B,  by 


and  set 


£=  1  1  , 
T)  -TJ 


H'=  B-'JfB, 


H'  =  B~XHB  , 


e  (z)  =  B-'6{z)  , 


Then,  equation  (1)  becomes 


k  (771)  =  B  1«c(«7i) 


H'$  (m)-H'e  (*,)=«  (*),  v  j7j  <  0 , 
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or,  in  component  form, 


A2  (^i)  —  Aj  0j  (t?i)  —  Kj  (»?l) 

Ai  e2  (T)l)  ~  A2  02  M  =*2  (Vi) 


V  Tji  <  0 


It  can  be  seen  from  above  analysis  that  H  and  H  can  be  diagonalized  simultaneously  by 
the  same  transform? ‘ion.  Therefore,  the  originally  coupled  equations  (1)  can  be  reduced 
to  the  uncoupled  equations  (8). 

If  we  express  the  ratio  A 3 /A2  as  having  the  following  dependence  on  0: 

Ai  _  1  -f  0 

A2  1  —  0 

then  the  parameter  0  must  be  expressed  as, 


y/h\2^2\ 

As  a  result,  the  solution  for  the  first  equation  in  (8)  can  be  obtained  as 

0i  (z)  _  J_  f  *1  (T)dT  1  \  rQx 

L{z)  2iriJc  A2£+(r)(r  -  2)  M 

O 

where  A  ( z )  is  an  arbitrary  entire  function.  C  is  a  contour  along  the  entire  branch  cut, 
and  extends  from  negative  infinity  to  the  interfacial  crack-tip.  The  function  L(z)  is  given 


where 


I(2)  =  2-i+*’+<t  , 


1  ,  1-0 

£  =  2?1"iTa’ 


and  ki  is  an  real  integer.  Integer  k\  is  chosen  so  that 


|  L(z)  |  =  0  (  |*|-  )  , 


as  |,zj  —»  0  , 


which  complies  with  the  restriction  of  equation  (2). 


Similarly,  we  can  obtain  that 


^  O 

02  (2) _ 1_  [  *2  (t)6.T 

T{z)  2 iri  Jc  Aj X+(r)(r  -  z) 


(11) 


v  . 

where  B  (2)  is  also  an  arbitrary  entire  function.  L  stands  for  the  complex  conjugate  of 


L. 


Returning  to  the  original  function  6(z), 


where 


and 


of  \  1  /  f  1  L(z)  1  L(z)  •  .  .1  dr 

0(z)  =  —  /  <  -rTTr^FK  r)  +  r  k(t)  ^ - 

4«  Jc  [ A2  L+(t)  A 1  Z,  (t)  j  r-; 


+  1(2)  A(2)C  +  I(2)  B(z)C 


r ,  1  ■> 

r  1  i 

1 

* 

1  — 

r  = 

V 

,  r= 

V 

v  1  . 

.  -  V  1  . 

C  =  ( i>  v  )T  >  C=  ( 1.  -i>  )T  • 


(12) 


III:  Some  asymptotic  results  of  the  Stieltjes  transform 

In  solving  the  Riemann-Hilbert  problem,  we  need  to  evaluate  the  integral 

W  =  f  •  (13) 

J—  00  fjl  —  2 

Setting  t  =  —rji,  we  get 

w = -  r  Jr**  ■  <i4) 

As  we  can  see  from  equation  (14),  —  I(z)  is  the  Stieltjes  transform  of  function  f(t). 
Here  we  want  to  study  the  asymptotic  behavior  of  the  Stieltjes  transform  as  z  — ►  0. 
Alternatively,  we  may  set  A  =  1/z  to  get 

I(z)  =  -\H[f;\)  ,  (15) 
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where 


os) 

Studying  the  asymptotic  behavior  of  (14)  as  z  — »  0  is  equivalent  to  studying  the  asymp¬ 
totic  behavior  of  (16)  as  A  — ►  oo. 


Suppose  that  /(f)  is  locally  integrable  in  (0,oo).  Recall  that  the  Mellin  transform  of 
/(f)  is  defined  by 

M(f;s]=  (17) 

JO 

and  set 

h(t)  = 


1  4*  t 


Then,  by  using  the  Parseval  formula,  we  can  obtain  that 

1  rr+ioo 

H[f;  A]  =  -T-.  /  A -M[h;  s]M[f\  1  -  s]ds  , 

Z7TZ  Jr—ioo 


(18) 


where  the  constant  r  is  such  that  Re(s)  =  r  lies  in  the  common  strip  of  analyticity  of  the 
Mellin  transforms  M[h;s]  and  M[f\  1  —  s). 


After  some  manipulations,  it  can  be  shown  that 


M[h-,s]  =  - 


sin  its 


(19) 


where  M[h;  s]  is  analytic  in  the  strip  0  <  Re(s)  <  1.  In  analogy  to  the  particular  problem 
of  interfacial  fracture  that  we  are  interested  in,  we  will  define  the  function  /(f)  as, 


/(«)  =  «“  (in  t)9  , 


(20) 


where  o-  =  ±2e,  or  0,  and  0  =  0,  or  1.  For  this  function,  the  Mellin  transform  M[/;s] 
only  exits  in  the  generalized  sense.  Let 

f(t),  t  €  (0, 1] 


m  = 


h(t)  =  { 


0,  t  €  (1,  oo) 


I  0,  t  €  (0, 1) 

/(f),  fe[i,oo) 
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Then,  we  may  write 


^;Al  =  l)(A)  =  f  mdi, 


>  =  1,2, 


Also  let 


Then, 


»[/iA)  =  i,(A)  +  i,(A). 


Gj(s)  =  M[h\ s]M[fj-,  1  —  s] , 


j  =  1,  2  . 


G(s)  =  1  -  s]  =  G,(s)  +  G,(s) 


In  addition,  from  the  Parseval  formula, 


L’w=hf:Zx"aMis' 


7  =  1,  2  , 


1 

W^]  =  oZ7E/  .  A-G^d 

Z7H  j  JTj~i oo 


Using  the  specific  function  /(£)  chosen  in  (20),  it  can  be  shown  that 

=  ~  [3  -  (l  +  ia)]W  sin  its  '  (26) 

In  the  above  we  can  see  that  Gi(s)  is  analytic  in  the  strip  0  <  Re(s)  <  1.  Since  M[/2;  1—  sj 
is  analytic  in  the  half  plane  Re(s)  >  1,  and  M [ft;  s]  can  be  analytically  continued  into  the 
entire  s-plane  as  a  meromorphic  function,  G2(s)  is  a  meromorphic  function  in  the  half 
plane  Re(s)  >  1  with  simple  poles  at  s  =  2,3,  *  *  *.  Then  in  equation  (24),  we  can  always 
choose  that  0  <  ri  <  1  and  r2  >  rj.  Observe  that  if  s  =  sx  +  ts2,  Gi(s)  has  the  property 


lim  Gi(si  +  is2)  =  0  , 

|j2|->oo 


<  Sj  <  r2  . 


Therefore,  we  can  apply  Cauchy’s  integral  theorem  to  equation  (25),  which  results  in 

# [/;  A]  =  £  res  {-A-Gt(s)}  +  f  A''G(s)ds  .  (2 

rj  <Re(*)<ra  2*1  Jr3-,oo 
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For  our  case,  it  is  easy  to  show  that  G(s)  =  0.  So  finally,  we  get 


«U;  A]  =  £  «* 

rj  <Re(»)<rj 


{- 


— } 

sin  its  J 


(1-1-  ia)]0+l  sin : 

Letting  r2  — ♦  +  oo,  we  get  an  infinite  asymptotic  series  for  H[f\  A]  as  A  — ►  oo, 

By  applying  the  above  analysis  to  our  particular  problem,  for  o  ^  0,  we 
following  asymptotic  results: 

?r2  cosh  7 ra  •  1  1 


l 


0  (-?.)“  In(-li) 


— oo  tjt  —  Z 


dr?i  = 


iir 


sinh  ira 


zta  In  z  — 


sinh2  to 


O* 


+0(  \z\  ) 


£ 


oo  1)1  —  Z 

ln(— 771) 


oo  T)i  —  Z 
0  i 


dr/! 

dJ/i 


+ 1  +  0(  M  ) 

smh  ttq  a 

5  (In  *f  +  j  +  0(  W) 


>  ,  as  z 


L 


oo  T)i  —  Z 


drji  =  Inz  -f  0(  \z\  ) 


(29) 


will  obtain 


-0.  (30) 
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Captions 


Figure  1:  Schematic  of  dynamic  growth  of  a  crack  along  a  bimateriai  interface. 

Figure  2:  Velocity  dependence  of  mismatch  parameter  r)  for  plane  stress  and  plane 
strain.  (Bimaterial  combination:  PMMA/steel) 

Figure  3:  Velocity  dependence  of  mismatch  parameter  0  for  plane  stress  and  plane 
strain.  (Bimaterial  combination:  PMMA/steel) 

Figure  4:  Velocity  dependence  of  mismatch  parameter  0  for  plane  stress  and  plane 
strain  at  the  vicinity  of  the  shear  wave  speed  of  PMMA.  (Bimaterial  combination: 
PMMA/steel) 

Figure  5:  Velocity  dependence  of  mismatch  parameter  c  for  plane  stress  and  plane 
strain.  (Bimaterial  combination:  PMMA/steel) 

Figure  6:  Velocity  dependence  of  the  real  part  of  mismatch  parameter  e  for  plane  stress 
and  plane  strain  at  the  vicinity  of  the  shear  wave  speed  of  PMMA.  (Bimaterial 
combination:  PMMA/steel) 

Figure  7:  Schematic  of  the  optical  set-up  for  CGS  in  transmission. 

Figure  8:  Selected  sequence  of  CGS  interferograms  of  a  growing  crack  in  a  three  point 
bend  interfacial  drop  weight  tower  experiment.  (Only  PMMA  side  of  PMMA/steel 
specimen  is  shown) 

Figure  9:  Velocity,  (a)  and  acceleration,  (b)  time  histories  for  the  experiment  shown  in 
Figure  8. 

Figure  10:  Time  histories  of  mismatch  parameter  e,  (a)  and  its  time  derivative,  (b)  for 
the  experiment  shown  in  Figure  8. 
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Figure  11:  Comparison  of  digitized  data  points  from  the  interferogram  corresponding 
to  t  =  9.5/zsec  in  Figure  8  with,  (a)  a  A'd-dominant  fit,  equation  (7.9);  (b)  a  higher 
order  transient  analysis  fit,  equation  (7.10).  (Crack  lies  along  the  negative  aq-axis) 

Figure  12:  Comparison  of  digitized  data  points  from  the  interferogram  corresponding 
to  t  =  23^sec  in  Figure  8  with,  (a)  a  /id-dominant  fit,  equation  (7.9);  (b)  a  higher 
order  transient  analysis  fit,  equation  (7.10).  (Crack  lies  along  the  negative  a^-axis) 

Figure  13:  Selected  sequence  of  CGS  interferograms  of  a  growing  crack  in  an  one  point 
bend  interfacial  gas  gun  experiment.  (Only  PMMA  side  of  PMMA/steel  specimen 
is  shown) 

Figure  14:  Velocity,  (a j  and  acceleration,  (b)  time  histories  for  the  experiment  shown 
in  Figure  13. 

Figure  15:  Time  histories  of  mismatch  parameter  e,  (a)  and  its  time  derivat.ve,  (b)  for 
the  experiment  shown  in  Figure  13. 

Figure  16:  Comparison  between  the  CGS  fringe  pattern  and  fitted  higher  order 
transient  stress  field,  equation  (7.10),  for  a  propagating  crack  in  a  PMMA/steel 
interface. 

Figure  17:  Velocity  time  history  for  the  experiment  shown  in  Figure  18. 

Figure  18:  Selected  sequer  ce  of  CGS  interferograms  of  a  growing  crack  in  an  one  point 
bend  interfacial  gas  gun  experiment.  (A  blunt  sta>  ter  notch  was  used) 

Figure  19:  CGS  interferograms  providing  visual  evidence  of  the  highb'  transient  nature 
of  dynamic  interfacial  crack  growth. 
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Material- 1 


Figure  1:  Schematic  of  dynamic  growth  of  a  crack  along  a  bimaterial  interface. 


( 1  ) 

Crack-Tip  Speed,  v/c  * 

Figure  2:  Velocity  dependence  of  mismatch  parameter  ij  for  plane  stress  and  plane  strain. 
(Bimaterial  combination:  PMMA/steel) 
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Figure  3:  Velocity  dependence  of  mismatch  parameter  j3  for  plane  stress  and  plane  strain. 
(Bimaterial  combination:  PMMA/steel) 
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Figure  4:  Velocity  dependence  of  mismatch  parameter  fi  for  plajie  stress  and  plane  strain  at 
the  vicinity  of  the  shear  wave  speed  of  PMMA.  (Bimaterial  combination:  PMMA/steel) 
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Figure  5:  Velocity  dependence  of  mismatch  parameter  e  for  plane  stress  and  plane  strain. 
(Bimaterial  combination:  PMMA/steel) 
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Figure  6:  Velocity  dependence  of  the  real  part  of  mismatch  parameter  c  for  plane  stress  and 
plane  strain  at  the  vicinity  of  the  shear  wave  speed  of  PMMA.  (Bimaterial  combination: 
PMMA/steel) 
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Figure  10:  Time  histories  of  mismatch  parameter  e,  (a)  and  its  time  derivative,  (b)  for  the 
experiment  shown  in  Figure  8. 
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Figure  11:  Comparison  of  digitized  data  points  from  the  interferogram  corresponding  to 
t  =  9.5/xsec  in  Figure  8  with,  (a)  a  /Codominant  fit,  equation  (7.9);  (b)  a  higher  order 
transient  analysis  fit,  equation  (7.10).  (Crack  lies  along  the  negative  .Ti-axis) 


Mixitv  79  degrees  ond  v/c 


Velocity  v/c 


Figure  12:  Comparison  of  digitized  data  points  fi  >m  the  interferogram  corresponding  to 
1,  -  2 3 //sec  in  Figure  8  with,  (a)  a  /\d-dominant  fit,  equation  (7.9);  (b)  a  higher  order  transient 
analysis  fit,  equation  (7.10).  (Crack  lies  along  the  negative  2q-axis) 
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Figure  14:  Velocity,  (a)  and  acceleration,  (b)  time  histories  for  the  experiment  shown  in 
Figure  13. 
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(a) 


(b) 

Figure  15:  Time  histories  of  mismatch  parameter  e,  (a)  and  its  time  derivative,  (b)  for  the 
experiment  shown  in  Figure  13. 
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Figure  [(*:  Comparison  between  the  CCS  fringe  pattern  ami  the  fitted  higher  order  transient 
stress  field,  equation  (7.10).  for  a  propagating  erark  in  a  PMMA/sleel  interface. 
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